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Abstract 


A rotorcraft roof sandwich panel has been redesigned to optimize sound power transmission loss (TL) 
and minimize structure-borne sound for frequencies between 1 and 4 kHz where gear meshing noise 
from the transmission has the most impact on speech intelligibility. The roof section, framed by a grid of 
ribs, was originally constructed of a single honeycomb core/composite facesheet panel. The original 
panel has coincidence frequencies near 700 Hz, leading to poor TL across the frequency range of 1 to 4 
kHz. To quiet the panel, the cross section was split into two thinner sandwich subpanels separated by 
an air gap. The air gap was sized to target the fundamental mass-spring-mass resonance of the double 
panel system to less than 500 Hz. The panels were designed to withstand structural loading from 
normal rotorcraft operation, as well as 'man-on-the-roof' static loads experienced during maintenance 
operations. Thin layers of VHB 9469 viscoelastomer from 3M were also included in the facesheet ply 
layups, increasing panel damping loss factors from about 0.01 to 0.05. Measurements in the NASA SALT 
facility show the optimized panel provides 6-11 dB of acoustic transmission loss improvement, and 6-15 
dB of structure-borne sound reduction at critical rotorcraft transmission tonal frequencies. Analytic 
panel TL theory simulates the measured performance quite well. Detailed finite element/boundary 
element modeling of the baseline panel simulates TL slightly more accurately, and also simulates 
structure-borne sound well. 
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1 Introduction 


Commercial rotorcraft are powered by drive systems comprised of complex transmissions, which 
contain sets of gears and shafts supported by bearings. As the gears rotate at high rates of speed, they 
induce vibration and noise in the transmission and throughout the rotorcraft. An example of a 
rotorcraft transmission region and roof section is shown in Figure 1. 

The goal of this work was to develop and evaluate acoustically tailored composite rotorcraft panels to 
reduce noise transmitted into the passenger cabin of a rotorcraft. The focus is on the structural roof 
panels, which are mechanically connected to the transmission, allowing strong gear meshing tones 
emanating from the transmission to pass into the panels and radiate into the cabin. Sound radiated by 
the transmission housing also impacts the ceiling panels acoustically, which transmit a portion of that 
sound into the interior. 

Composite materials are sometimes used to construct lightweight stiff panels for rotorcraft, which 
reduce weight, but also lead to increased sound radiation into the rotorcraft due to their reduced 
impedances and increased sound radiation efficiencies. Trim panels constructed of layers of foams and 
thin plates are often attached to the panels to reduce sound transmission, but are expensive, bulky, and 
heavy. The scope of this work is to design composite fuselage panels which do not require trim panels, 
and actually reduce, rather than increase, the noise from transmission tones. 



Figure 1. View of the inside of the transmission region in a commercial rotorcraft. 

A sample of in-flight test data which illustrates key frequencies of interest for noise reduction is shown 
in Figure 2. While several tones are present, two frequencies were identified based on sound quality 
assessments as most critical in this example - those of the main rotor bull gear mesh (at around 1060 
Hz) and the input pinion gear mesh (at around 3100 Hz). Modeling and measurements of the baseline 
and optimized panels will therefore focus on frequencies between 1 and 4 kHz. 
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Figure 2. Sample of in-flight sound pressure data for a commercial rotorcraft. 


In this work, we have developed and demonstrated technologies and methodologies for designing 
composite fuselage panels which radiate less sound into the rotorcraft interior. We constructed and 
tested two 4-ft by 4-ft roof panels- one using current composite manufacturing approaches, and a 
second panel based on the optimized acoustic technologies applied under this project. The program has 
achieved the following specific overall objectives: 


1. A baseline, 4-ft by 4-ft, flat, composite, beam-stiffened panel representative of the roof panel in a 
commercial helicopter was designed. The baseline panel meets the same requirements, including 
weight and structural integrity, as the roof panels in commercial helicopters. 


2. The baseline composite panel was fabricated using standard manufacturing techniques used at Bell 
Helicopter. 

3. Methods for designing composite panels were developed that reduce the vibrational energy 
propagating into and through the panels while improving the sound radiation characteristics. 

4. An optimized composite roof panel was designed and fabricated to minimize vibration transmissibility 
and sound radiation over frequencies spanning the dominant gear meshing tones of a typical rotorcraft 
(1-4 kHz). Practical constraints (weight and structural integrity) were included in the optimization. 

This work provides NASA and the rotorcraft community with design methods that can be used to 
acoustically tailor composite panels for rotorcraft. This work also provides NASA with representative 
composite panels that can be used to evaluate future treatment concepts. 
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2 Baseline Composite Panel 


2.1 Design and Construction 

A baseline composite panel with geometry and properties representative of those of a Bell Helicopter 
429 roof structure has been designed and constructed. A cross-section within the desired region of the 
429 roof assembly was selected as the basis for the panel design, as shown in Figure 3. Equivalent 
composite materials were chosen to replace the existing metallic materials for the upper and lower 
skins, edge band laminate, and core. Thicknesses and surface densities of the materials in the center 
sandwich panel are listed in Table 1, and equivalent structural material properties are discussed in 
Section 2.3. Note that Table 1 includes a layer of viscoelastomer attached to the inner surface which is 
typical of that bonded to many Bell roof panels to mitigate noise transmission. 

The upper and lower skins of the center sandwich panel are made of plies of Cytec G30-500/5276-1 
Carbon/Epoxy Fabric. A layup of plies was used which provides the best equivalent stiffness, as 
measured by the product of Young's Modulus and Moment of Inertia, to that of the 429 core-stiffened 
metallic panel. The core thickness is the same as that of the production panel, so that equivalent 
stiffness is controlled by the panel face sheet properties. Face sheet stiffness is defined using: 

Em^m — EqAc and G M A M — G^A^ , 

where E is the extensional modulus (nominally 8.2 Msi for the Cytec fabric), G is the shear modulus 
(nominally 400 ksi), A is the area, M indicates metallic properties and C indicates composite properties. 
Since the area is the product of width and thickness (/), and the panel maintains the same width, these 
equations may be rewritten as: 

Em^m = E c t c and Gm^m = ■ 

Several composite laminates were evaluated to find the layup with the best match to that of the 
production panel. A three ply layup with orientation angles (in degrees) of [0/45/0] was eventually 
chosen to provide the best combination of extensional and shear stiffness in comparison to the metallic 
panel. 

A similar analysis was performed for the edge band laminate, except that instead of shear or extensional 
stiffness, equivalent bending stiffness is the critical parameter. The equation for determining equivalent 
bending is: 

= E c t c 3 . 

The layup that best satisfies this constraint is [0/45/45/0/0]s, where 'S' indicates symmetry about the 
center ply for the remaining angles. The edge, therefore, is constructed of 10 total layers of Cytec. 

The last component in the composite baseline panel design is an equivalent composite core material. It 
is not practical to achieve a stiffness comparable to that of an Aluminum core with a composite core, so 
the next critical driver for design is strength. Hexcel Kevlar core with a density of 47 kg/m 3 (2.9 lb/ft 3 ) 
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and a mean Shear modulus of 96.5 MPa (14 ksi) was found to have the best equivalent strength to that 
of the metallic core. 

Design loads for the core stiffened panel area are critical for the upper skin in compression. Upper skin 
applied ultimate loads are set to 150% of the flight limit loads, and include components in both in-plane 
directions, and in-plane shear. Design loads for the edge laminates are also based on flight load limits, 
and include in-plane forces as well as a moment applied about the fore-aft direction. A margin of safety 
greater than 1 is required for both core and edge fabric laminates based on an elevated temperature 
wet open hole compression allowable. Other design constraints are (a) local skin wrinkling stability, (b) 
core crush in the ramp area, and (c) core shear due to a man-load on the roof panel. The baseline panel 
meets all of these requirements. 

Standard size aluminum I-beams were chosen with cross-sectional properties that best matched the 
properties of the variable cross-section production roof beams and transverse intercostal beams. In the 
baseline panel, the roof beams are 1016 mm (40 inches) long, and the intercostal beams are 762 mm (30 
inches) long. The flange widths of the roof and intercostal beams are 76.2 mm (3 inches) and 50.8 mm 
(2 inches) respectively. Both beams are 102 mm (4 inches) high, with 3.96 mm (0.156 inch) flange and 
web thicknesses. The roof beams and transverse intercostal beam webs are connected by aluminum 
shear clips, and the top flanges are connected by four aluminum splice straps that provide an 
economical representation of the joints at the transmission mounting points. The beams, shear straps 
and splices are shown in Figure 4. The panel design is detailed in Figure 5 and Figure 6. The roof beams 
are connected to the panel with 5/16" diameter titanium protruding shear head pins and titanium 
collars spaced at a nominal 6.5 Diameter pitch. These fasteners and spacing are representative of the 
majority of the Bell model 429 roof beam to panel fasteners. 

Figure 7 shows the composite portion of the baseline panel construction, and Figure 8 shows the fully 
assembled panel. The total panel mass is 15.8 kg (34.7 lb), with 5.64 kg (12.4 lb) in the composite, 10.1 
kg (22.3 lb) in the beams, and 0.60 kg (1.33 lb) in the fasteners and shear straps. The completed panel is 
coated with an Epoxy VOC Compliant Primer (MIL-P-85582, TY I) Light Green (water-base). 
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Table 1. Cross-sectional materials, thicknesses, and surface densities for center region of baseline panel. 


Material 

Thickness 

(in) 

Thickness 

(cm) 

Surface 

density 

(lb/in 2 ) 

Surface 

density 

(kg/m 2 ) 

Paint 



0.00006 

0.042 

Face sheet 

0.0079 

0.0201 

0.000443 

0.311 

Face sheet 

0.0079 

0.0201 

0.000443 

0.311 

Face sheet 

0.0079 

0.0201 

0.000443 

0.311 

Adhesive 


0.0000 

0.00035 

0.245 

Honeycomb 

0.500 

1.2700 

0.00048 

0.671 

Adhesive 


0.0000 

0.00035 

0.245 

Face sheet 

0.0079 

0.0201 

0.000443 

0.311 

Face sheet 

0.0079 

0.0201 

0.000443 

0.311 

Face sheet 

0.0079 

0.0201 

0.000443 

0.311 

Paint 



0.00006 

0.042 

Visco 

0.04 

0.1016 

0.00222 

1.56 

Total 

0.5S7 

1.49 

0.00618 

4.67 




Figure 3. Cross-section of Bell 429 Roof Assembly, dimensions in inches. 
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Figure 4. Baseline panel viewed from above (top) and below (bottom). 
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Figure 5. Stackup of skins and honeycomb core. 



Figure 6. Baseline panel geometry (not to scale). 
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Figure 7. Composite portion of baseline panel. 
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Figure 8. Fully assembled baseline panel viewed from above (top) and below (bottom). 
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2.2 Structural Material Properties 

A finite element (FE) model was constructed of the baseline panel to analyze its sound transmission 
characteristics. An original FE model had several errors, including inaccurate and incomplete material 
properties, incorrect orientation of the core ribbon and warp directions, and incorrect beam 
thicknesses. The errors led to the FE model overestimating resonance frequencies by 20-30%. To better 
determine the actual structural material properties, two 4" x 8" coupon panels cut from a 12" x 12" 
newly constructed panel intended to be identical to that of the baseline panel were tested statically to 
compare static stiffnesses, as well as dynamically to compare measured and simulated resonance 
frequencies of the fundamental panel modes. The panels were cut so that the length of one was aligned 
along the ribbon direction of the core, and the length of the other aligned with the warp direction. 

Since the test panels are quite small and extremely light, traditional modal testing with attached 
accelerometers was impractical, since the mass loading from the accelerometers would have been 
significant. Instead, the panels were suspended from wire adjacent to a loudspeaker, which ensonified 
the surfaces with airborne white noise. A Laser Doppler vibrometer (LDV) measured the resulting panel 
surface vibration, shown in Figure 9. Four modes of vibration are visible for each panel - two near 2 
kHz, and two near 3.5 kHz. 

FE models of the test coupons were generated, and the modal frequencies were compared to those 
extracted from the measurements. The model was initially corrected to include the mass density of the 
layers of adhesive (Cytec FM 300K film adhesive, .05 I b/ ft 2 , estimated to be nominally .008 inches thick, 
http://www.cemselectorguide.com/pdf/FM 300 081211.pdf) , as well as the paint. The stiffnesses of 
the paint and adhesive are ignored. The adhesive is assumed to rigidly connect the face sheets and 
core. Adjusting the surface mass density led to improved agreement between measured and simulated 
resonance frequencies, but the FE model was still about 10-15% too stiff. This meant that either or both 
of the face sheet and core stiffnesses provided by the material suppliers were incorrect. 

Bell conducted static measurements of the core and face sheet materials, and updated the material 
properties. Notable changes include: 

o a heavier core, with a 10% increase in density; 

o a softer core, with about a 20-30% decrease in stiffness; and 

o softer face sheets, with about a 10% reduction in stiffness. Also, Bell discovered that 
the laminate stiffness differs slightly (about 6%) for the 'tool' and 'bag' sides. 

The updated material properties greatly improves the agreement between the measured and simulated 
test coupon modes (within +-5%). The final structural material properties for the laminates and core are 
listed in Table 2. The compression and tension moduli, as well as the tool and bag side moduli, were 
averaged to compute the values used in structural FE dynamic analyses. The variation in the 
tool/bag/compression/tension moduli is small - about 4%. Also, the in-plane Poisson's ratio for the face 
sheet is a smeared estimate based on the stack of three plies. A key lesson learned from this experience 
is that FE modeling and test coupons are invaluable for insuring material properties are properly 
characterized. 
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Table 2. Final laminate and core material properties used for analysis. Directions 1 and 2 are in-plane, and 3 is 
through the thickness. Note the face sheet Poisson's ratio is smeared over the full set of three plies. The edge 
material has material properties identical to those of the face sheets. 


Table 2a - Selected center panel material properties 


Property 

Face Sheets 

Kevlar Core 

En, E 22 (GPa/Msi) 

57/8.3 

- 

V12 

0.21 

- 

Gi3 (MPa/ksi) 

- 

139/ 20.1 (ribbon) 

G 23 (MPa/ksi) 

- 

68 / 9.8 (warp) 

yo(kg/m 3 /lb/in 3 ) 

1550/0.0560 

47/0.0017 


Table 2b - 

Edge panel properties 

Property 

Value 

En, E 22 (GPa/Msi) 

54/7.8 

V12 

0.21 

yo(kg/m 3 / lb/in 3 ) 

1550/0.0560 

t (mm / in) 

2.0/. 079 



Frequency (Hz) 

Figure 9. Vibration of test coupons ensonified with acoustic field. The length of panel 1 is aligned with the 
ribbon direction of the core, and the length of panel 2 aligned with the warp direction of the core. The LDV 

normal velocity units are uncalibrated. 
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2.3 Structural Modeling (Finite Elements) 

2.3.1 Bell Structural FE model 

A structural FE model of the baseline panel was constructed at Bell, as shown in Figure 10. All 
components are modeled with quadratic solid elements, with a 0.5" spacing along the width and length. 
The upper and lower skin plies were each discretely modeled with a single element through the 
thickness. The core was modeled with one element through the thickness. The edge band laminate was 
modeled with six elements through the thickness. A cross-section of the panel is shown in Figure 11. 

The I-beams were modeled with a single element through the thickness, and are connected to the panel 
in the locations of the fasteners using CBUSH (spring) elements to represent the bolts. The stiffnesses of 
the spring elements are based on the bolt material and cross-section. The beam cross-section and 
CBUSH elements (shown as red circles) can be seen in Figure 12. 

Initially, the model was created with elements with an aspect ratio which conforms to Bell standard 
practices. However, this led to an unacceptably large number of degrees of freedom, so the model was 
then remeshed with a 20:1 aspect ratio through the thickness (the element edge length along the width 
and length of the panel is 20 times that of the element thickness). The transition area from the upper 
and lower skins to the edge band region was simulated by discretely modeling all plies in the edge band 
and then tapering the plies to the upper and lower skins as shown in Figure 13. 



Figure 10. Finite element model of baseline panel. 
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Upper Skin 


Core 


Edge Band 



Figure 11. Panel cross-section near edges. 
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Figure 13: Final panel cross-section near edges. 

2.3.2 Penn State Structural-Acoustic Model 

The Bell FE model resolves all face sheet plies individually, and is suitable for structural integrity analysis. 
Penn State analyzed Bell's FE model of the baseline panel using NX/NASTRAN. Unfortunately, the model 
required over a day to compute normal modes for frequencies up to 1 kHz. Since the project requires 
analysis of frequencies up to 4 kHz, and the Bell model cannot be analyzed easily for those frequencies, 
Penn State generated a lower-resolution FE model for acoustics analysis. The lower resolution model 
represents the face sheets with fewer through-thickness elements, reducing the model size by over 60%. 
Quadratic solid elements are still used, however, consistent with Bell's model. Here are the details of 
the modeling reductions: 

• The edge laminate regions (without honeycomb core) are now modeled with two through- 
thickness elements, with aggregate material properties computed based on Bell's orientations 
and material properties (see Table 2). 

• The inner panel face sheets are modeled with a single element on the inner and outer surfaces, 
rather than three, with aggregate material properties computed based on Bell's orientations 
and material properties (see Table 2). 

• The honeycomb core is modeled with a single through-thickness element. 

Based on Penn State's experiences with similar structures, little error was anticipated with using 
aggregate face sheet properties rather than modeling each ply layer individually. The Aluminum beams 
are modeled with the same mesh resolution as Bell's model, and are connected to the panel using the 
same spring definitions, representing the bolts. 

In addition to Bell's bolts, Penn State also applied a spatial stiffness in the transverse and in-plane 
directions to represent the frictional coupling of the beams and panel, which are compressed together 
by the bolts. The interface stiffness is based on an assumed percentage of the faces of the mating 
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materials faces in contact, which is related to an assumed surface roughness. We use guidance from 
Garvey [1], who measured the effects of preload on the overall elastic moduli of stacked motor laminate 
sheets. Our baseline structure differs from stacked common laminate sheets, however, since the beam 
and panel materials are different. We therefore use the stiffnesses of the softest of the two materials in 
contact - the composite, as this will be the limiting value. Since the surface roughness is not known, we 
assume the interface elastic moduli are 5% of those of the composite moduli (1.3 Msi in compression 
and 0.41 Msi in shear) and assume an interface thickness of 0.079 in. To convert to an array of springs 
between the nodes of the beams and the panel, the stiffnesses are combined with the ratios of the 
interface areas and the number of springs between interface nodes. 

2.4 Acoustic Modeling (Boundary Elements) 

To compute the sound radiated by the panel and its effects on panel radiation damping, an acoustic 
boundary element (BE) model was constructed. Every structural element has a corresponding acoustic 
element. The overall model is assumed to be baffled by the adjacent fuselage sections (or by walls in a 
Transmission Loss Test Chamber). The fluid loading and sound radiation were computed using the 
lumped parameter approach of Koopmann and Fahnline [2]. The resulting fluid loading matrices were 
applied to the FE model using the component mode synthesis approach in ARL/Penn State's CFIAMP 
(Combined Flydroacoustic Modeling Programs) analysis procedure [3, 4], 

2.5 Measurements 

Several types of vibroacoustic measurements were conducted on the baseline panel, including 
experimental modal analysis, surface averaged vibration response functions, radiated sound power, and 
sound power transmission loss. The modal analyses were conducted for free and clamped boundary 
conditions. The free boundary condition data were used to validate and guide updates to the FE model. 
The clamped data were used to confirm the edge boundary conditions applied to the panel when 
mounted in the NASA SALT facility [5]. Finally, the sound power and transmission loss (TL) 
measurements made in SALT validate the overall FE/BE model, and act as a baseline for future optimized 
panel performance assessments. 

See Appendix A for detailed test procedures for these measurements. 

2.6 Vibro-acoustic behavior 

2.6.1 Effective Flexural Wavespeeds 

The effective flexural wavespeeds are examined to better understand the nature of sound transmission 
through the panel. A cross section of a honeycomb core, composite facesheet panel is shown in Figure 
14. The wavespeeds of honeycomb core/composite facesheet sandwich panels are dominated by 
moment effects at low frequencies, where the facesheets contract and expand in flexure, and by shear 
effects at high frequencies, where the core rigidity resists transverse motion. The effective bending 
wavespeed Cb may be determined from: 
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2N 


( 1 ) 


c 2 =- 


pM \(Ps h ) 


Ap s hN 2 
orD 


where D = + , N = G core h( 1 + 1 / hf , and G core = yjG^G^ . a is the overall panel mass 

density, h is the honeycomb core thickness, iV is the shear rigidity, D is the flexural rigidity, E is the 
Young's Modulus, t is the overall facesheet section thickness, vis Poisson's ratio (in the in-plane, or '1-2' 
direction spanning the overall facesheet sections), and G are the transverse shear moduli in the 
different directions of the honeycomb 'weave' (ribbon and warp). Note the importance of core 
thickness to D, and correspondingly to q>. 

The critical frequency, where the effective bending wavespeeds match the sound speed in air (c 0 ), is: 


co. 


<P,h 1 

d i _ c ipy 

N 


( 2 ) 


Using the material properties in Table 1 and Table 2, the baseline panel flexural rigidity D is nominally 
3200 N-m, the mean shear rigidity Af varies from 0.94 to 1.93 MN-m (the shear modulus is different in 
the warp and ribbon directions), the facesheet section effective Poisson's ratio is 0.21, and the overall 
mass density per unit area is 3.11 kg/m 2 (this includes the mass of adhesive and paint, but not of any 
added external layers of viscoelastomer). Figure 15 compares the center and edge panel bending 
wavespeeds computed using these properties, with the edge panel (no honeycomb core) in-plane 
Poissons ratio of 0.25. The mean shear and effective bending wavespeeds, along with upper and lower 
bounds based on the variable shear moduli in the in-plane directions, are shown. The center panel 
bending waves become supersonic at frequencies above about 700 Hz, with the edge panel section 
waves remaining subsonic up to 5 kHz. Also, the effective bending wave speed begins to approach the 
shear wave speed upper limits above 5 kHz. 



Figure 14. Schematic of typical honeycomb core/composite facesheet panel. 
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Figure 15. Analytic wavespeeds for the center panel and edge panels; the mean and upper and lower bounds on 
the shear and effective sandwich bending wavespeeds are shown. 

2.6.2 Resonances 

The reduced version of the FE model does not explicitly model every composite ply, but groups them 
into combined layers to reduce model size. To ensure the simplified version is accurate, Figure 16 
compares a set of low frequency modes extracted from both Bell's original model and Penn State's 
model with free boundary conditions. The modal frequencies are nearly identical. The modes of the 
panel may be subdivided into categories - central panel (with stiffening honeycomb core), as shown in 
Figure 17, and edge panels (composite material only), as shown in Figure 18. There are also modes of 
the beam structure (Figure 19). The low-order center panel mode shapes resemble simply supported 
panel modes, with the ribs acting as the simple supports. Also, the outer edge mode shapes resemble 
those of clamped-free panels. Figure 20 compares the measured and simulated center panel resonance 
frequencies, and shows that the FE model captures the frequencies within 10% uncertainty. 

The resonance frequencies and panel dimensions may be used to estimate flexural wavespeeds in the 
center panel by assuming edge boundary conditions. For example, the wavespeed c m „ of a given mode 
shape at its resonance frequency C0m n may be inferred from: 


where k m „ depends on the assumed boundary conditions. The center panel is connected to the beam 
stiffeners via a transition region between the honeycomb core and pure facesheet material, as shown in 
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Figure 21. To assess whether the boundary condition best emulates simply supported or clamped 
conditions, the modal wavespeeds were computed assuming: 



+ 


f nn ^ 1 


u J 


(4) 


for simply supported boundary conditions and 


k =, 


(2m + l)^Y f(2n + l)^ 2 


2a 


+ 


2b 


(5) 


for clamped conditions, where a is the width (nominally 38 inches) and b is the length (nominally 33 
inches). The flange half-widths were added to the center panel width to compute the nominal a and b 
dimensions for the wavenumber analysis since the flange is thin and moves with the center panel. 

Figure 22 compares inferred modal wavespeeds to those computed using analytic formulae for infinite 
panels, and shows that the panel modes are effectively simply supported. The figure also shows the 
range of flexural wavespeeds in the panel due to the differing honeycomb core stiffnesses in the ribbon 
(stiffest) and warp (weakest) directions. The data are shown again in Figure 23, but are limited to 
modes with only a half wave across either the length (m=l) or width (n=l). The m=l modes feature 
bending waves which travel predominantly along the width, or ribbon (stiffer) direction, and the n = 1 
modes are aligned primarily along the length, or warp (weaker) direction. Figure 23 shows that these 
modes follow the analytic curves fairly well, particularly the n = 1 modes. 

This exercise shows that the center panel region is indeed effectively simply supported, with minor 
mode shape length and width variations. This means that simple analytic formulas for simply supported 
panels may be used to guide optimized panel design. 



Figure 16. Comparison of low-frequency modal frequencies extracted from Bell rigorous FE model and Penn 
State approximate FE model; diagonal line indicates perfect agreement (slope of 1.0). 
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( 1 , 1 ) 

127 Hz FE/122 Hz Exp 


( 24 ) 

251 Hz FE/247 Hz Exp 



( 34 ) ( 2 , 2 ) 

486 Hz FE / 487 Hz Exp 426 Hz FE / 438 Hz Exp 

Figure 17. Low-order center panel modes of baseline panel, FE simulations and FE and experimental resonance 
frequencies. Values in parenthesis indicate the number of half-waves along the length and width of the panel 

respectively. 



Figure 18. Example of mode shapes in edge material; left image zooms on edge material. 
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Beam in-plane bending, 593 Hz 


Beam flange twisting, 668 Hz 


Figure 19. Examples of beam modes of baseline panel, FE simulations. Comparable experimental mode shapes 
were not identified since the beams were not part of the measurements. 



Figure 20. 


Measured vs. simulated center panel resonance frequencies for baseline panel, 'n' corresponds to 
the mode order along the width direction of the center panel. 
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Figure 21. Connections between beam and center and edge composite panels. 
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Figure 22. Flexural wavespeeds of the center panel inferred from resonance frequencies and assumed modal 
wavenumbers for simply supported and clamped conditions. The inferred wavespeeds are compared to analytic 

estimates in the warp and ribbon honeycomb directions. 
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Figure 23. Same as Figure 22, except that only modes with a half wave in the length (m=l) or width (n=l) 

directions are shown. 


2.6.3 Damping Loss Factors 

Overall damping loss factors were extracted from the experimental modal analysis data, and are shown 
in Figure 24. The center panel modes are labeled in the figure, and are highly damped near 500 Hz, 
slightly below the acoustic coincidence frequency. Above about 2 kHz, all modes have damping loss 
factors between nominally 0.006 and 0.035. Figure 25 shows the center panel damping only, along with 
estimates of the damping induced by sound radiation using the relation: 


7 had 


~ ° rad 


P„ C o ' 

copy 


( 6 ) 


where cw is the sound power radiation efficiency of the center panel modes, p 0 is the mass density of 
air, and c 0 is the sound speed in air. At and above coincidence, we assume a ra d is 1 (as shown in the 
acoustic BE calculation in Figure 26), and see that the estimated radiation loss factor is generally an 
upper bound of the center panel mode loss factors. However, the actual center panel modes, 
particularly at higher frequencies, couple strongly with the beams and edge panels, as seen in Figure 17 
and Figure 19. This coupling increases the effective modal masses of those modes. Increasing the panel 
mass density by 50% to approximate this effect reduces the radiation damping, and provides a better 
match with the higher frequency panel modes above 800 Hz. Rigorous calculations of individual modal 
masses would provide improved agreement, but this simple comparison is sufficient to illustrate the 
importance of radiation damping on the center panel vibration. It is critical, therefore, that acoustic BE 
(or similar) modeling be included in any panel vibration calculation to account for radiation damping. 
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Figure 24. Measured modal damping loss factors. Center panel modes are distinguished by their (m,n) mode 

orders 
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Figure 25. Measured center panel modal damping loss factors vs estimates of radiation damping. 


24 




5 



Frequency (Hz) 

Figure 26. BE calculation of radiation efficiency of baseline panel for diffuse drive in dB (0 dB indicates radiation 

efficiency of 1). 

2.6.4 Forced Response Simulations and Measurements 

Two forced response calculations were used to better understand the critical sound transmission paths 
through the baseline panel and to guide design optimization efforts: a structural drive at one of the 
beam joints (mimicking transmission drives as shown in Figure 27), and a diffuse acoustic field applied 
over the surface of the center panel (shown in Figure 28), as well as both the center and edge paneling. 
All radiated sound transfer functions were computed using ARL/Penn State's CHAMP analysis capability 
[3, 4] (see Figure 29). Based on the loss factor measurements, we assume a global mechanical loss 
factor of 0.01, and apply the acoustic BE model to the structural FE model to capture the effects of 
radiation damping. 

The simulated and measured point mobilities for three center panel drives are averaged and compared 
in Figure 30. The measurements and simulations agree well, and are also mostly bounded by the infinite 
panel mobilities above 1 kHz. Upper and lower bound infinite panel mobilities are computed using: 


G mf = 


CO 

%pM 


(7) 


and the variable core shear moduli in the ribbon and warp directions. The upward trending of the 
conductance is caused by the increasing dominance of the core shear modulus with increasing 
frequency. The mobilities above 1 kHz (the important frequency range for this project) are influenced 
by many modes, with only modest variability. This high modal overlap, along with all modes being 
supersonic (with unit radiation efficiency) indicates that any quieting strategy must target nearly all 
modes equally, rather than a few strong sound radiators. 
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The sound power radiated by the panel due to excitation by a diffuse acoustic loading field is measured 
in NASA Langley's SALT [5] facility, and simulated using the CHAMP software. To simulate the boundary 
conditions in the SALT facility, the outer edges of the panel are fixed in the FE model. Only the panel 
regions (center and edge) are driven with the diffuse loading; the beams are unloaded. Figure 31 shows 
the simulated sound power radiated by the center panel and the edge region for frequencies above 400 
Hz (the lower frequency limit of the NASA reverberant room), and shows that, as expected, the faster 
bending waves in the center panel are responsible for much higher transmitted sound than the subsonic 
waves along the edges. 

The powers in Figure 31 are compared to the usual assumption of incident sound power in a reverberant 
room: 


P =■ 

me 


C n W- 


( 8 ) 


where w,„ is the reverberant room energy density and S is the panel surface area. Using the blocked 
pressure assumption we can approximate w in as: 


w- 



(9) 


where p is the acoustic pressure at the boundary. Since we apply a unit pressure loading to the panel, 
the squared pressure must be unity. 

The transmission coefficient is then 


r rad 


S/4p 0 c 0 


( 10 ) 


where P m d is computed by CHAMP. The virtual transmission loss (VTL) is the inverse of the transmission 
coefficient expressed in dB: 


VTL = 10 log 


r S/4p 0 c 0 ' 


r rad 


(ii) 


The computed VTL for the center panel region is compared in Figure 32 to that of the analytic infinite 
panel estimate in Fahy and Gardonio [6], integrated over all angles of incidence: 


k/2 
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More details on the VTL procedure are in Appendix B. 
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The FE/BE-based VTL is similar to the analytic estimate, establishing confidence in the modeling. 
However, the small TL of the panel shows that sound sources are poorly attenuated, and that panel 
redesign to mitigate transmission noise is clearly warranted. Finally, Figure 33 shows that the overall 
VTL (which includes the center and edge panels) is for most frequencies within 3 dB of that measured in 
the NASA SALT facility, establishing further confidence in the modeling procedures. The coincidence dip 
of the edge paneling near 5.5 kHz is not as strong in the measured data, likely because the edge paneling 
is not a large continuous section, but a narrow frame, such that simplified infinite panel TL techniques 
are no longer accurate near coincidence. 

Figure 31 and Figure 33 also include analytic estimates of radiated sound power and TL from infinite 
panel representations of both the center and edge regions. The total panel analytic TL estimates are 
computed by summing the total power radiated by the total areas of each region. The analytic 
estimates show that the edge paneling radiates more power above about 4 kHz. Optimizing the TL of 
the center panel, therefore, will only significantly improve overall panel TL below 4 kHz. 

Finally, radiated sound power simulations and measurements for a transverse structural drive at a rib 
joint are compared in Figure 34. The measurements were made using an intensity probe on the 
anechoic side of the SALT facility. As with the diffuse drive, the simulated and measured data are 
generally within 3 dB (except at 500 Hz and 1 kHz, which show up to 5 dB differences). 



Figure 27. Typical rotorcraft installation. Roof panel is in green below the strut. Yellow arrow denotes 

oscillatory transverse force. 


27 



Figure 28. Diffuse acoustic loading applied to center region of baseline panel. 



Figure 29. ARL/Penn State's CHAMP analysis capability [3]. 
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Figure 30. Simulated and measured drive point mobilities, averaged over three drive locations. Infinite panel 
upper and lower bound mobilities computed for core ribbon and warp directions. 
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Figure 31. Simulated input and radiated sound power. 
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Figure 32. TL for center panel only, FE/BE VTL vs. analytic simulations. 
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ure 33. Measured vs. simulated TL for center and edge regions for diffuse acoustic drive. 
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Figure 34. Measured vs. simulated sound power for structural drive at beam joint. 
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3 Panel Optimized for Reduced Noise 


3.1 Design Optimization Goals 

Design options were pursued which balance vibroacoustic reductions with weight and structural 
integrity constraints. The center region of the panel was adjusted to minimize structure-borne and air- 
borne sound transmission without adding excessive weight. In this context, net weight is defined as the 
offset between additional weight of the optimized panel (as compared to the baseline panel) and the 
weight of soundproofing insulation and/or damping materials that can be eliminated as a result of 
increased vibro-acoustic performance of the optimized panel: 

Net Weight = Optimized Panel Weight - Baseline Panel Weight - Insulation/Damping Eliminated Weight 

Typical surface weights range from 1 - 2.4 kg/m 2 (0.2 - 0.5 lb/ft 2 ) for standard soundproofing packages 
and 3.4 - 9.7 kg/m 2 (0.7 - 2 lb/ft 2 ) for VIP applications. One objective of the design optimization was to 
reduce or eliminate soundproofing materials based on the vibro-acoustic performance of the optimized 
panel. The following generic soundproofing configuration would typically be applied to the baseline 
panel to improve vibro-acoustic performance: 

1. 60% coverage constrained layer damping material, 1.56 kg/m 2 (0.32 lb/ft 2 ) 3M Aearo 
Technologies C-2206-03PSA, 1 mm (0.04 in) thick, and 

2. 2.54 cm (1 in) thick, 6.7 kg/m 3 (0.42 lb/ft 3 ) Johns Manville Microlite AA fiberglass insulation 
blanket. 

During the optimization study, the cross-sectional materials used in the center panel were varied, while 
ensuring that the flexural rigidity was maintained, and that material strength constraints were met (see 
Section 2.1 for a discussion of structural design constraints). An important design constraint is to ensure 
symmetric laminates with a minimum of 10% fibers in each direction. Also, lamina data factored with 
open hole 'knockdown factors' conservatively approximates Bell's typical allowable design curves. 
Typical critical design allowable is for an open or filled hole in compression at an elevated temperature 
and moisture soaked condition. 

3.2 Quieting Technologies 

Based on the analysis results to date, the baseline panel represents a challenging noise reduction 
problem. The center panel bending waves are supersonic at frequencies above about 700 Hz, so that all 
modes of vibration radiate sound strongly throughout the key frequency range of interest (1-4 kHz). 

There are several design changes that could reduce sound transmission of the panel, including: 

attaching treatments to the panel; 

reducing the shear stiffness of the core material to reduce the effective panel flexural 
wavespeeds, increasing the coincidence frequencies and reducing sound power radiation 
efficiencies; 

adding structural damping to the panel, reducing vibration and therefore radiated sound power; 
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introducing structural barriers adjacent to the roof beams to reduce the structure-borne energy 
flow into the center panel; and 

splitting the panel into two subpanels separated by an air gap. 

The design changes must not (a) add significant weight or thickness, and (b) reduce structural integrity. 

3.2.1 Traditional Treatments Attached to Panel 

There are currently two typical methods of reducing cabin noise - applying thin sheets of constrained 
layer damping (CLD) material and/or layer(s) of fiberglass insulation to the interior face of the panel. 

The damping material reduces vibration amplitudes at panel resonances, thereby reducing radiated 
sound. The mass of the insulation reduces the sound power transmitted through the center panel, 
particularly at high frequencies. 

Common off-the-shelf CLD material is ISODAMP C-2206 from E-A-R Aearo Technologies (owned by 3M). 
The material is 1 mm thick with a surface density of 1.56 kg/m 2 , and according to the vendor, provides 
an added damping loss factor of 0.115 to a thin sheet of Aluminum at 1 kHz. The added mass to the 
panel by the CLD is significant, since the bare panel surface density is only 2.48 kg/m 2 . This added mass 
will also slow down the panel flexural waves, but not to the point where the coincidence frequency is 
increased enough to reduce radiation efficiency above 1 kHz. 

A commonly used acoustic blanket is the Microlite AA from Johns Manville. The lightest Microlite AA 
blanket weighs 0.168 kg/m 2 (small compared to the panel surface density), and provides about 4.5 dB of 
sound attenuation at 1 kHz and about 12 dB at 3 kHz. The blankets are mounted within a thin plastic 
covering which is attached to the edges of the panel. 

3.2.2 Reduced Stiffness Core Materials 

The classical approach for increasing sound power transmission loss in honeycomb sandwich panels is to 
slow down the shear waves in the core material to the point where they are subsonic [7, 8]. The Kevlar 
core used in the baseline panel has an effective mean shear wave speed of 821 m/s, which is well above 
the air sound speed of 344 m/s. Reducing the mean core shear modulus from 97 MPa to 15 MPa would 
reduce the shear wave speed to 287 m/s, eliminating the panel coincidence frequency such that all 
flexural/shear waves are subsonic. Using infinite panel theory, the TL curves in Figure 35 show the 
potential effects of slowing down the core shear waves. 

The results show that the upper bound for sound power transmission loss is dictated by simple mass law 
theory. Further reductions of airborne sound power are only possible with added noise barriers, such as 
the acoustic blankets commonly used by industry, or by adding other structural panels, with an 
accompanying airgap, normal to the panel. Also, it is unlikely that a panel with a core with a mean shear 
modulus of 15 MPa would meet the structural integrity constraints described in Section 2.1. 

In an alternative approach investigated by NASA Langley Research Center, voids may be machined into 
the core material to reduce the transmission of energy through the panel [9]. This approach, however, 
reduces the panel overall bending stiffness, and may not be acceptable for a rotorcraft roof panel. 
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Figure 35. TL for baseline panel, and optimized panel with soft core and subsonic wave speeds. 

3.2.3 Embedded Viscoelastic Treatments 

Thin layers of viscoelastomer may be embedded within the facesheet sections, or within the core itself. 
Peters et. al. [8] reported loss factors of 0.05 and greater with core-embedded viscoelastomer, and 
showed a 3-5 dB transmission loss improvement over more lightly damped panels. Reducing radiated 
sound at low frequencies via structural damping improvements will be difficult due to the dominance of 
the radiation damping (in cases of high core shear stiffness). To investigate the potential benefits of 
added structural damping, the structural loss factor was increased from 0.01 (measured in the baseline 
panel) to 0.05 over all frequencies in the FE model. Figure 36 and Figure 37 show simulated reductions 
in radiated sound power for structural and acoustic diffuse drives on the FE/BE model. As expected, 
sound power reductions due to added structural damping are most pronounced at higher frequencies, 
with a mean reduction of 2-3 dB, and a reduction at resonance peaks of 4-6 dB. These results are 
promising, and consistent with those of Peters [8]. 
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Figure 36. Effects of added structural damping on radiated sound power due to structural drive at beam 

intersection. 



Figure 37. Effects of added structural damping on radiated sound power due to acoustic diffuse drive. 
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3.2.3.1 VHB Viscoelastomer from 3M 


The FE/BE results indicate that inserting layers of thin viscoelastic material in the face sheets, as well as 
between the beam stiffeners and panel, will lead to improved acoustic performance. The face sheets 
used in helicopter sandwich panels, however, are extremely thin and lightweight (each ply in the 
baseline panel is 0.0079 inches thick). Replacing one of the layers, in this case the center ply, with a 
layer of viscoelastomer requires similarly thin and light damping material. 3M corporation 
recommended using their VHB 9469 adhesive, which is formulated to have high damping properties 
near room temperature and at frequencies between 1 and 4 kHz. Replacing the center layer of each 
face sheet with a layer of the VHB material (which is nominally 0.005" thick) leaves outer and inner 
carbon layers with 0 and 90 degree ply orientations. There are no 45 degree ply orientations in the new 
design, however. Structural integrity calculations described in Section 3.3.1 confirm the acceptability of 
this approach. 

Figure 38 shows the temperature- and frequency-dependent shear modulus and loss factor for the VHB 
material. Shear Moduli (computed assuming a Poisson's ratio of 0.499 and Young's Moduli from DMA 
testing) and loss factors are compared at 20 and 30 degrees C in Figure 39. The loss factors are quite 
high, ranging between 0.7 and 1.1 between 1 and 4 kHz. Though promising, the damping benefits of the 
VHB material were checked by performing experimental modal analyses on two test coupons 
constructed by Bell. 
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Figure 38: Shear modulus and loss factor nomograph for 3M VHB 9469 (provided by 3M). 
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Figure 39: Young's modulus and loss factor for 3M VHB 9469 at 20 C and 30 C. 

3 . 2 . 3. 2 Coupon Testing and Analysis 

Two test coupons constructed with layers of the 3M VHB 9469 material sandwiched between sheets of 
carbon fiber plies (the same material used in the baseline panel) were tested using experimental modal 
analysis at Penn State. The panel dimensions (19" x 23") were chosen to avoid modal degeneracy, so 
that each structural mode is distinct in frequency and easily identified. The same Hexcel Kevlar 0.5" 
thick core (1/8" cell size, 3.3 lb/ft 3 ) used in the baseline panel was used for the test coupons. The two 
coupons were constructed using different approaches. The first panel has facesheets with cocured 
carbon fiber and VHB, such that part of the VHB fused with epoxy in the carbon fiber sheets. This 
formed a hybrid structure with uncertain properties. A second panel with pre-cured carbon fiber sheets 
post-bonded with the VHB was therefore also constructed. A non-destructive inspection (NDI) of the 
post-bond panel was conducted to ensure complete adhesion. An image of the NDI test is shown in 
Figure 40, and indicates excellent contact between the plies and VHB, with no voids (blue colored 
regions) detected. 
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Figure 40: NDI scan of improved post-bonded face-sheet/VHB/face-sheet layup. Small blue squares are pads 
used for testing, and not voids. White regions are above sensing range, but still indicate excellent adhesion. 

Complex modes were extracted from experimental modal analysis data, and loss factors and resonance 
frequencies were compared for the two panels. Figure 41 shows a few of the measured mode shapes, 
which are clean and symmetric, providing confidence in the quality of the panel construction. Figure 42 
compares the modal loss factors for the two panels for frequencies up to about 5 kHz. The post-bond 
approach consistently yields higher damping, and both construction approaches lead to strong damping 
improvements at 1 and 3 kHz, where the dominant transmission tones occur. Based on these data, the 
post-bond approach was used for the optimized panel, and a nominal panel structural damping value of 
0.05 is assumed for acoustic performance simulations. Flowever, this construction technique is costlier 
and more time consuming than a co-cured approach. Bell will weigh performance and cost should it 
decide to pursue the optimized panel concept in future rotorcraft. 

Replacing the center carbon layer reduces the face sheet stiffness, thereby reducing the flexural 
wavespeeds. The measured mode shapes may be used to estimate modal wavenumbers, which 
combined with modal frequencies determine modal wavespeeds (see Appendix C for details on the 
approach). Figure 43 shows the inferred modal wavespeeds for the two panels along with a least- 
squares data fit using thin plate theory and the wavespeed of the baseline sandwich panel. The baseline 
panel properties were used to infer the modulus for the optimized panel facesheets, which have the 
center ply replaced with 3M VFIB 9469 adhesive. Using the baseline sandwich panel properties, it was 
found that the sample panel wavespeed reasonably matches the modal wavespeed estimates if the 
sample panel facesheet wavespeed is about 80% of the baseline panel facesheet wavespeed and the 
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shear wavespeed is the same. To infer the sample panel facesheet modulus, thin plate theory was 
used along with the ratio of the wavespeeds: 


CjJ^sample _ a/A/PA _ V^A 2 / 1 ~ ^ )P, _ Q q (13) 

^bjs, baseline \f^b 1 P h /l2(l — U^)p t 

The panel thicknesses and Poisson's ratio are approximately the same while the ratio of the densities is 
about 0.75. The final estimate of the facesheet modulus then becomes E s =0.3Eb. a significant 
reduction. The reduced stiffness and wavespeed will lead to higher acoustic coincidence frequencies. 
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Figure 41: Selected measured mode shapes and loss factors of post-bond test coupon. 
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Figure 42: Measured structural loss factors for two test coupons with embedded VHB 9469. 
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The panels also provide an opportunity to verify the FE modeling procedure for sandwich panels with 
layers of VHB material. Figure 44 shows a schematic of the cross-sectional modeling of the panels. The 
coupon was modeled with 4,370 parabolic solid elements. Each ply layer was modeled with one 
element through its thickness, and four elements represent the Kevlar core. The ribbon direction was 
modeled along the length of the panel. The adhesive layers between the inner plies and the core were 
not modeled explicitly, but the layer masses were simulated instead by increasing the adjacent ply mass 
densities from 0.311 kg/m 2 (0.000443 lb/in 2 ) to 0.558 kg/m 2 (0.000793 lb/in 2 ). The final modeled (0.757 
kg) and measured (0.753 kg) masses match almost exactly. 

The viscoelastomer Young's Moduli were varied over several center frequencies per the data shown in 
Figure 39. The Poisson's ratio was set to 0.499. Modes of the panel with free boundary conditions were 
extracted using NASTRAN for each property set, and modal frequencies were determined based on 
proximity to the center frequency of each set. Figure 45 compares the measured and simulated 
resonance frequencies, which agree to within +-4%. Figure 46 compares measured and simulated 
structural loss factors, which agree well for frequencies above 1 kHz. Below 1 kHz, the simulated loss 
factors are higher than the measured ones. Flowever, since this project focuses on frequencies between 
1 and 4 kHz, we have not pursued the cause of this discrepancy. Overall, the good agreement between 
measured and predicted resonance frequencies and loss factors confirm the modeling procedure and 
the underlying material properties. 
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Figure 44: FEM model for test coupons with embedded VHB material. 
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Figure 45: FE vs. measured resonance frequencies for test coupon with embedded VHB 9469 and free boundary 

conditions. 
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Figure 46: FE vs. measured damping loss factors for test coupon with embedded VHB 9469. 


3.2.3.3 Projected performance improvement with embedded VHB 9469 

The structural damping of the baseline panel FE/BE model was increased from 0.01 to 0.05 based on the 
test coupon measurements. The acoustic sound power transmission loss calculations were repeated, 
and compared to the baseline panel TL in Figure 47. The increased damping leads to about a 4 dB noise 
reduction, which while helpful, is not spectacular. Other means of improving noise transmission are 
therefore still required. 
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Figure 47: CHAMP simulations of the effects of VHB on baseline panel. 


3.2.4 Band-Gap/Phononic Crystal Structure-borne Sound Barriers 

Phononic crystals, or arrays of structural discontinuities, can mitigate structure-borne sound 
transmission of waves in certain frequency ranges. A formulation for the exact solution for multiple 
scattering of flexural waves in a thin plate has been obtained. A technical report that contains all the 
details of the formulation is provided in Appendix E. The problem consists of an array of inclusions, of 
an arbitrary number and arrangement, embedded in a thin plate, and the multiple scattering of a planar 
incident wave is analyzed. The key assumptions in the formulation include: the plate is thin (hence the 
classical thin-plate theory is used), uniform in its mechanical and geometrical properties, and has an 
infinite extent. Three types of inclusions (scatterers) have been included in the formulation: 

• Voided: mathematically this assumes that the perimeter of the inclusion is free of forces and 
moments. Physically this represents drilled through-thickness holes in the plate. 

• Rigid: mathematically this assumes that the perimeter of the inclusion is clamped, and hence 
both the displacement and slope vanish. 

• Elastic: this assumes that inclusions are thin plates. Mathematically this requires the 
continuation of the displacement and slopes across the interface. Although the classical thin- 
plate theory may not be valid for the inclusions when the diameter of the inclusions is in the 
same order as their thickness, it is assumed that the classical thin-plate theory remains valid as 
the first order approximation. 
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The formulation has been implemented as computer code so that simulations for the wave transmission 
effects with various material and geometric configurations can be performed. In the simulations here, 
the inclusions are considered rigid. Some preliminary simulation results are presented and discussed in 
the following. In this set of simulations, only a simple geometry is considered: the scatterers are all 
identical, of a radius of 0.015 m (30 mm diameter), and are spaced 0.05 m apart, in a square 
arrangement. The simulations are run at a few discrete frequencies of interest: 1 kHz, 3 kHz, 6 kHz, and 
10 kHz. 

3.2.4.1 Voided Inclusions 

The amplitude of the wave field in the vicinity of the voids is shown in Figure 48. In this figure, the upper 
left corner is for/= 1 kHz, the upper right corner is for/= 3 kHz, the lower left corner is for/= 6 kHz, and 
the lower right corner is for/= 10 kHz. The color scale has a maximum of 2.5 (red), and minimum of 0 
(blue). 



Figure 48: Amplitude of acoustic pressure surrounding the void inclusions at four frequencies: 1 kHz (upper-left), 
3 kHz (upper-right), 6 kHz (lower left), and 10 kHz (lower-right). 

It can be seen from this figure that the scattering effects are extremely weak at all four frequencies. 
There are only slight traces of reflection at 3 kHz and 6 kHz. 
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3. 2.4.2 Rigid Inclusions 

Similarly, Figure 49 shows the amplitude near the rigid inclusions at the same four frequencies. In this 
case the scattering effects are very strong. In fact they are so strong that wave transmission is 
prohibited almost entirely. And they are effective in all the four frequencies computed. In other words, 
there is virtually no frequency dependency. In essence, the scatterers effectively form a wall that blocks 
the transmission of the wave. 

3.2.4.3 Elastic Inclusions 

For elastic materials, there are wide ranges of values for different material properties. A series of 
simulations are performed in which the inclusion's flexural rigidity and mass density change 
proportionally relative to those of the plate, which represent a most common characteristic in natural 
materials. Results of two most representative cases are presented here. In the first case, the inclusion's 
flexural rigidity is 100 times of that of the plate, and the mass density is 10 times of the plate. This case 
is designated as "hard elastic" case. In the second case, the inclusion's flexural rigidity is 0.01 times of 
that of the plate, and the mass density is 0.1 times of the plate. This case is designated as the "soft 
elastic" case. The resulting wave fields at those four frequencies near these elastic inclusions are shown 
in Figure 50 and Figure 51. It can be seen that neither case provides an effective barrier for the wave 
transmission. Furthermore, it appears that a hard elastic inclusion is not as strong a scatterer as a soft 
elastic inclusion. 



Figure 49: Amplitude of acoustic pressure surrounding the rigid inclusions at four frequencies: 1 kHz (upper-left), 
3 kHz (upper-right), 6 kHz (lower left), and 10 kHz (lower-right). 
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3. 2.4.4 Discussion of Simulation Results 

This set of cursory (preliminary, rather than systematic) simulations indicates that, except the case of 
rigid inclusions, the scattering effects in all other cases are very weak. Aside from the obvious reason 
that the inclusions in those cases are "weak scatterers," the more important reason is that the 
arrangements do not match the conditions that are needed for band gap formation. Typically a 
phononic band gap forms when the wavelength in the host (the thin plate) approaches the spacing 
between the adjacent scatterers. In the base sandwich panel, at 3 kHz, the wave speed is around 650 
m/s (according to Figure 15), which gives a wavelength of ~0.22 m. This wavelength is much greater 
than the spacing and hence in the low-frequency regime, the inclusions provide very little impedance. 
Even at 10 kHz, the highest frequency in the simulations, the wavelength is ~0.14 m, still far greater than 
the spacing. 



Figure 50: Amplitude of acoustic pressure surrounding the soft elastic inclusions at four frequencies: 1 kHz 
(upper-left), 3 kHz (upper-right), 6 kHz (lower left), and 10 kHz (lower-right). 
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Figure 51: Amplitude of acoustic pressure surrounding the hard elastic inclusions at four frequencies: 1 kHz 
(upper-left), 3 kHz (upper-right), 6 kHz (lower left), and 10 kHz (lower-right). 

There are two directions that can be explored to address this issue and to bring a band gap into the 
picture. One is to reduce the sound speed in the sandwich panel. The potential challenge with this 
direction is that the changes thus incurred to the sandwich panel compromise the structural stiffness 
requirements. The second direction is to enlarge the spacing. For example, for targeting the noise 
reduction at 3 kHz, the scatterers can be arranged with a spacing of ~0. 25-0. 27m. The potential 
problems with this direction are that an increasing the spacing usually would correspondingly require a 
larger inclusion size; and that the overall size of the sandwich panel (0.9m X 0.9m) may not allow a 
sufficient number of inclusion for the band gap to form. Of course, it is also possible to explore the 
combination of these two directions. 

3.2.4.5 Proposed Design to Rigidity the Inclusions 

Probably the most feasible, and if successful, the most effective, way to introduce the band gap effect 
for suppressing the wave transmission is to take a closer look at the "exception case": the rigid 
inclusions. 

In many theoretical studies, an elastic solid can be used to approach void or a rigid body by setting its 
material properties to extreme. However, in the set of simulations that has been performed for elastic 
inclusions, setting the flexural rigidity of the elastic inclusion to an extremely high vale (such as 10 9 times 
of the thin plate) does not bring a scattering effect that even remotely resembles the rigid inclusion 
case. 
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A closer look reveals the following underlying reason: the boundary conditions for the rigid inclusions 
are a mathematical idealization, and actually are not a real physical rigid inclusion. For a physical 
inclusion of extremely high rigidity but having the same thickness of the plate, the inclusion may 
undergo rigid body translations and rotations, collectively called rigid body motions. However, the 
mathematical boundary conditions for the so-called rigid inclusion prevent such rigid-body motions 
completely. This set of mathematical boundary conditions can be physically realized as having the 
inclusions not only rigid, but also fixed in space. This insight prompted the proposal of the following 
mechanism to prohibit or impede rigid body motions of the inclusions. 

In the proposed design, called stiffened inclusions, the inclusions are much longer than the thickness of 
the plate so their ends protrude out of the thin plate. These ends are bound to a stiffening panel, one on 
each end of an inclusion. Finally, to reduce the added weight, the center portions of the stiffening panels 
are eliminated. A 3D rendering of this conceptual design of such stiffened inclusions is shown in Figure 
52. 



Figure 52: 3-D solid model rendering of "stiffened inclusions" concept design. 

In this proposed structure, because the ends of the inclusions are all tied together, the rigid body 
rotation of individual inclusions is essentially eliminated, and the only possible motion is the rigid body 
translation. The rigid body translations of individual inclusions are further greatly impeded because 
those stiffening panels, along with the inclusions, form another sandwich panel with significantly higher 
rigidity (because of much larger thickness), and in essence, eliminates the rigid body translation. 

3.2.4.6 Band-gap target frequency range 

The stiffened inclusion concept was evaluated using the baseline panel design. There are two main 
frequencies of interest near critical transmission noise tones - 1 and 3 kHz. Typical mode shapes of the 
baseline panel near those frequencies are shown in Figure 53. While the structural wavelengths are too 
long for a practical inclusion array to be effective, a 3 kHz array seems feasible. To minimize interaction 
of the array with nearfield evanescent waves near the beams, an offset of a few inches is required. 
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Acoustic wavelengths are ~4/5 that of structural wavelengthsat 1 kHz 



Acousticwavelcngthsare 7/16 that of structural wavelengthsat 3 kHz 


Figure 53: Structural and acoustic wavelength estimates at 1 and 3 kHz. Yellow rectangle indicates potential 

region for stiffened inclusion array. 

3.2.4.7 Effective Panel Properties 

The multiple scattering analysis/simulation is based on the classical thin-plate model, in which the thin 
plate is assumed to be uniform and isotropic. However, the sandwich roof panel behaves more like a 
thick panel. Fortunately, we can still use the plate scattering analysis simulation procedure near a target 
frequency (3 kHz for this application) by computing an equivalent flexural rigidity at that frequency. The 
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parameters that characterize the wave propagation are the thickness h, the mass density p and the 
flexural rigidity D of the plate. The total thickness and the mass density are calculated based on the 
physical configuration of the panel, giving h = 13.9 mm and p = 223.3 kg/m 3 . The equivalent flexural 
rigidity is calculated according to sandwich panel bending wave speed c/, based on the following relation 
in the thin-plate model: 


D = 



(14) 


where co= 2nf, and /is the wave frequency. The thus-retrieved flexural rigidity is shown in Figure 54. The 
sonic crystal will be targeted to have a stop band centered around 3000 Hz. From Figure 54, a flexural 
rigidity of 836 N-m is chosen for the thin-plate model, which gives a wave speed of 556 m/s at 3000 Hz. 
The wave speed in the thin-plate model for frequencies up to 8000 Hz is plotted in Figure 55, where the 
wave speed based on thick-plate model is also plotted for comparison, showing a good match at 3000 
Hz. 



Figure 54: Equivalent flexural rigidity from the wave speed based on a thick-plate model. 
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Figure 55: Wave speed for the thin-plate model based on the equivalent flexural rigidity at 3000 Hz (red), 
compared with estimate based on the thick-plate model (blue). 

3.2.4.8 Initial Design 

Commercial off-the-shelf inserts, such as those shown in Figure 56, may be used as the array elements. 
Such inserts provide the entire surface of the through hole in the inserts as an adhesion surface. The 
weaknesses of this choice are the added weight and the significant reduction of the radius of the rods, 
reducing the rigidity of the overall inclusion assemblage. In numerical simulations, the inclusion's radius 
will be the radius of the insert, which is somewhere between those of the flanges and the housing. In 
this choice, the inclusion will be solid steel or fiber reinforced composite rods. 



Figure 56: Various forms of potted inserts for sandwich panels. (Taken from the web site of the provider: 
http : //www . clipnuts . com/potted_in_inserts . html) 
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Some of the parameters used in the simulation are selected in advance: 


• At least two rows of inclusions are required based on preliminary scoping studies. Three rows can 
achieve a better performance but takes up significantly more of the panel's surface. 

• The spacing of the inclusions is between 50.8 and 152.4 mm (2 and 6 inches). Denser packing 
would enhance the wave blockage effect, but risks reducing the overall structural integrity of the 
panel. Larger spacing would reduce the wave blockade effect. 

• The radius of the inclusions is between 5 and 15 mm (diameter: from 10 to 30 mm). 

Two initial designs, both of square arrangements, are assessed for the optimized panel 

• Design A: r = 15 mm, d = 56 mm; and 

• Design B: r = 15 mm, d = 100 mm. 

where r is the radius of the rigid inclusion, and d is the side length in a square unit cell. The difference 
between the two designs is that, physically, Design A targets the first band gap (3 kHz); whereas Design 
B targets the second band gap (1 kHz). Their transmission spectra for the baseline panel are shown in 
Figure 57. The band gap central frequency for Design A (red) is already near 3000 Hz (actually 3500 Hz). 

Next, the arrangements are finely tuned such at the central frequency is located closer to 3000 Hz. 



Figure 57: The forward displacement power ratios for the two initial designs. Red: Design A; Green: Design B. 
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3. 2.4.9 Tuning Design A 


The radius of the inclusion is varied while the spacing remains unchanged. The radius r is varied in both 
directions, ranging from 13 mm to 17 mm in 1 mm increments. The resulting spectra are shown in Figure 
58. The variation of the gap's central frequency with the radius is shown in Figure 59. It is also noted 
that as r increases, the bottom of the gap appears very flat. 



Figure 58: Forward displacement spectra with different inclusion radius: 
blue: r = 13 mm, red: r = 14 mm; black: r = 15 mm; purple: r = 16 mm; turquoise: r = 17 mm. 



Figure 59: Variation of the central frequency of the band gap with the inclusion radius r. 
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Similarly, while r remains unchanged (at r = 15 mm), the spacing d is varied, from 56 mm to 64 mm in 2 
mm intervals. The resulting spectra are shown in Figure 60; and the variation of the central frequencies 
with the spacing d is shown in Figure 61. 



Figure 60: Forward displacement spectra at different spacing at a fixed inclusion radius r = 15 mm. Blue: d = 56 
mm; Red: d = 58 mm; Black: d = 60 mm; Purple: d = 62 mm; Turquoise: d = 64 mm. 



Figure 61: Variation of the band's central frequency with spacing while the inclusion radius is fixed at r = 15 mm. 
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It can be seen that adjusting eithe r or d can fine tune the band gap. In general, a larger incluson radius 
results in a lower transmission coefficient. Thus, r = 15 mm is selected, and the spacing, based on a 
simple interpolation in Figure 61, is chosen as d = 63 mm. 

3.2.4.10 Tuning Design B 

The radius of the inclusions is varied, from 15 mm to 20 mm in 1 mm increments, while the spacing is 
fixed at d= 100 mm. The resulting spectra are collected in Figure 62. The variation of the central 
frequency with the radius is shown in Figure 63. It can be seen that shift of the central freqency is rather 
limited by varying the radius. Since the gap has a central frequency near 4700 Hz, it woud require a 
significant decrease in radius to shift the band gap to the left. This would deteriate the performance of 
the band gap. In other words, it is impractical to move the band gap by adjusting the radius alone. 



Figure 62: Forward displacement spectra with different inclusion radius: black: r = 15 mm, red: r = 16 mm; 
green: r = 17 mm; blue: r=18 mm; purple: r=19 mm; turquoise: r = 20 mm. 
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Figure 63: Variation of the central frequency of the band gap with the inclusion radius r. 

Increasing the spacing would shift the central frequency to the lower range. Since this is the only 
parameter to be adjusted, variable step sizes are used. The calculated spacings are d = 100 mm, 110 
mm, 115 mm, 120 mm. By 120 mm, the central frequency is below 3000 Hz. The collected spectra are 
show in Figure 64. The variation of the central frequency with spacing is shown in Figure 65. 



Figure 64: Forward displacement spectra with different spacing while the inclusion radius is kept at r = 15mm. 
Black: d = 100 mm, red: d = 110mm; green: d = 115 mm; blue: d =120 mm. 
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Figure 65: Variation of the central frequency of the band gap with the inclusion spacing d. 

Using spacing d = 120 mm, the central frequency is around 2900 Hz. According to Figure 64, increasing 
the radius by 1 mm, the central frequency increases approximately by 200 Hz. Thus, the final parameters 
for design B are chosen as: r = 15.5 mm; d = 120 mm. 

3.2.4.11 Final Designs 

The forward displacement spectra for the two final designs are compared in Figure 66. In plotting these 
curves, a finer frequency step size is used. Based on the spectra alone, both designs offer significant 
wave blockage. At the lower point, the transmission coefficients are approximate 0.45 and 0.25, 
respectively. This is the transmission coefficient for the displacement. In terms of power, they would be 
0.21 and 0.06; or, 6 and 12 dB reductions. 
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Figure 66: Forward displacement spectra for the final designs. Red: Design A. Green: Design B. 

The panel has an overall size of 762 mm x 941 mm (30"x36") between the beams. Leaving a space of 100 
mm between the loading beam and the sonic crystal, this leave 562 mm x 741 mm that can be occupied 
by the sonic crystal. For Design B, this would leave 5 inclusions in the short side and 7 on the long side. 
The second row would have 3 inclusions in short side and 5 inclusions on the long side. This essentially 
occupies the entire area. This design is not practical. 

For Design A, there are 10 inclusions in the shorter side, and 12 inclusions in the longer side. The second 
row would have 8 inclusions in the shorter side and 10 inclusions on the longer side. This seems to be a 
reasonable design. The overall layout is shown in Figure 67. 
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Figure 67: Overall layout of the sonic crystal (Design A) on the 30"x36" area of the panel. The black border 
represents the 30"x36" (762 mm x 914.2 mm) area. The distance between the sonic crystal and the shorter side 
is 110.6 mm, and the longer side is 97.5 mm. The spacing is 63 mm. 

3.2.4.12 Evaluation of Commercial Insert Geometry 

A metallic insert was identified, Shur-Lok Part #5169-S-8, as shown in Figure 68, with a shank outer 
diameter of 20.6 mm, and enlarged ends of diameter of 23.0 mm to allow wrenching holes; and a flange 
diameter of 41.3 mm. Using a slightly larger diameter is to account for the increase of the effective 
radius due to the flange, a diameter of 23.0 mm is used in the simulations. All other modeling 
parameters remain the same. 
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Figure 68: Diagram and geometry of the selected insert from Shur-Lok catalog (size 8 is used here). Diameter £ is 

used in the simulations. 

Since one of the key variables, the radius of the inclusions, has been fixed, simulations are run when the 
other primary design parameter, the spacing of, is varied from 40 mm to 60 mm in 5mm increments. The 
resulting spectra for the displacement amplitude in the forward direction are show in Figure 69. It can 
be seen that the central frequency of the band gap occurs at around 3400 Hz for d = 55 mm (brown 
curve); and around 2600 Hz for of = 60 mm (dark purple). When the spacing is small, the first band gap is 
located at a much higher frequency. 
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Figure 69: Spectra for different scatterer spacing: red: d = 40 mm, green: d- 45 mm; blue: d = 50 mm; light blue: 

d = 55 mm; dark purple: d = 60 mm. 

With interpolation, it is expected that a spacing of 57 mm would be about right. The spectrum is 
calculated in a finer step size in the same frequency range for d = 57 mm and is shown in Figure 70. The 
numerical data show that the central frequency of the band gap is located between 2980 and 3000 Hz in 
a computation using 20 Hz as the step size. 



Figure 70: Spectrum for the final design: r = 11.5 mm, d =57 mm. 
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3.2.4.13 Practical Design Considerations 

The retaining panel concept shown in Figure 52 is not practical for a commercial rotorcraft roof panel 
due to weight and space constraints. Alternate approaches include attaching the retaining frame 
directly to the outside surfaces of the inner and outer face sheets, effectively creating a 'doubler' 
picture-framed region over the insert array. Also, the honeycomb core in the insert region may be 
realigned so that the cells are perpendicular, rather than parallel to, the inserts. The doublers and the 
realigned honeycomb core both serve to mitigate rotational motion near the inserts. It is doubtful, 
however, that a perfectly irrotational array can be practically constructed. 

3.2.4.14 Effects of Finite Boundaries on Wave Scattering 

The effects of finite boundaries on wave scattering by an array of inclusions may be approximated by 
using additional wave sources to represent reflections from boundaries. A point source is produced by a 
concentrated time harmonic lateral force acting at a particular point on the plate. A series of such 
forces acting along a line segment forms a line source, which, by incorporating different phases, can 
represent a clamped edge or free edge. Furthermore, a series of not so closely located point sources 
could also mimic the effect of a row of rivets that holds a stiffening rib to a panel. 

Comparing Point Sources with Planar Incident Wave 

The first set of simulations explores the feasibility of using point source as a wave source, and also helps 
validate the sonic crystal design. For this, a side-by-side comparison is made for an array of two rows of 
inclusions. On one side, a planar incident wave is the wave source; on the other; a series of point 
sources are lined up parallel to the rows of the inclusions. Figure 71 shows an example of such 
comparison, in which the displacement amplitude is shown at the target frequency of 3000 Hz. The left 
figure has a planar wave as the source; and the right figure has an array of 17 equally spaced point 
sources, of the same strength and phase as the source. The point sources are spread out along a line 
segment that is noticeably longer than the inclusion array. This ensures an almost-planar wave form 
impinging onto the sonic crystal, which allows the comparison with the planar incident wave case. Note 
that each of the point source results in a singular wave field. In the computations, the amplitude of the 
incident wave due to these wave sources is normalized by the amplitude of the source 10 mm away. 

In a series of simulations, the distance from the point source array to the sonic crystal is varied. The 
simulations show that the wave field passing through the sonic crystal is essentially unchanged. The 
distance affects the reflected wave field in the space between the point source array and the sonic 
crystal, largely due to the change in the pattern for the standing wave between the line source array and 
the inclusion array. In the case of the planar incident wave, there is no standing wave but just 
interference of the two waves, the incident and the reflected wave, traveling in opposite directions. This 
observation confirms that the filtering effects of the band gap due to the inclusions are independent of 
the incident wave forms. 
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Note that due to the singular nature of the point loads, the normalization schemes used in the two cases 
are different and hence the direct comparison of values (or the exact colors) between the two cases are 
not meaningful. All comparisons focus on the pattern and the relative displacement amplitude within 
each figure. 



Figure 71. Distribution of displacement amplitude in vicinity of rigid inclusions at frequency /=3000 Hz. Left: 
planar incident wave. Right: incident due to an array of point sources 

Two Sonic Crystal Arrays 

In this set of simulations, two sonic crystals, each consisting of two rows of inclusions in square 
arrangement, are embedded in the panel, and one array of point sources is placed on the outside of 
each sonic crystal. In essence, the configuration adds a mirrored image of the sonic crystal and point 
source arrays as in the right image in Figure 71. 

In the first of this set of simulations, the location of the point source arrays is varied. Figure 72 
compares the displacement field of the entire panel for different locations of the source array, identified 
by the distance to the center of the panel. The sonic crystals are at fixed locations, with the sonic crystal 
rows located 271.5 mm and 328.5 mm away from the center of the panel, but with the point source 
arrays located 350 mm, 400 mm, 450 mm and 500 mm, respectively, from the center. 

Figure 73 shows the distributions of the displacement amplitude along the x-axis (horizontal symmetry 
line of the panel) for the four cases shown in Figure 72. In this figure, the locations of the inclusions are 
marked by the magenta vertical dashed lines; and the locations of the point source arrays are marked by 
the vertical dashed lines in the color matching that of the curve. The wave fields inside the two sonic 
crystals are essentially the same, with the exception of the first case in which the source arrays are 
located extremely close to the sonic crystal. Again the result suggests that the source location is an 
unimportant factor; and that the sonic crystals exhibit desired filtering effects at the target frequency of 
3000 Hz. Due to the space limit on the panel, there is not enough room to fit in a third row of sonic 
crystals. Otherwise, the filtering effect could be improved. 
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Figure 72. Displacement amplitude in the panel at 3000 Hz for the four locations of the point source arrays. 
Top left: 350 mm; top right: 400 mm; bottom left: 450 mm; bottom right: 500 mm. 
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Figure 73. Distribution of relative displacement amplitude along x-axis for the cases shown in Figure 72. 
Vertical dashed lines represent the locations of the source arrays (in matching color), and inclusions 
(magenta). Legend: red: 350 mm; green: 400 mm; blue: 450 mm; and orange: 500 mm. 
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Full-Panel Simulation 


To observe the behavior of the panel with the embedded sonic crystal, in the next series of simulations, 
the full panel is used. The sonic crystal now consists of two rings of inclusion. Based on the design of the 
sonic crystal, and the available space in the panel, four configurations as shown in Figure 74 are 

simulated. In the full panel, the load is transmitted into the panel through a set of I-beams that form a 

rectangular enclosure of dimensions 762 mm x 914.4 mm (30 in x 36 in). The four configurations have 
11, 12, 13 and 14 inclusions respectively in the horizontal directions, and 11 inclusions in the vertical 
direction, on the outer ring. The inner rings have two less inclusions in each direction. Two arrays of 
point sources are located 914.4 mm apart, one on each side of the sonic crystal. Within each array, 17 

point sources are equally spaced 50 mm apart, aligned in parallel with the sonic crystal. 
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Figure 74. Layouts of four configurations. In each configuration, the outline represents the portion of the roof 
panel within the I-beams through which the load is transmitted into the panel. Top left: llxll inclusions in the 
outer ring. Top right: 12x11 inclusions in the outer ring. Bottom left: 13x11 inclusions in the outer ring; Bottom 

right: 14x11 inclusions in the outer ring. 

Figure 75 shows the displacement field over the entire panel when the point source arrays are excited at 
the target frequency of 3000 Hz. Similarly, Figure 76 shows the distribution of the displacement 
amplitude along the x-axis in the panel. This set of simulations show a very different scenario compared 
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to Figure 72 and Figure 73. The most prominent difference is the strong resonance pattern formed 
inside the enclosure of the sonic crystal. This suggests, since the filtering cannot entirely block the wave 
from transmitting into the space enclosed by the sonic crystal, that the transmitted wave can resonate 
under appropriate conditions (when wavelengths align with the span of the geometry). 



Figure 75. Displacement amplitude in the panel at frequency 3000 Hz for the four configurations shown in Figure 

74. 
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Figure 76. Distribution of relative displacement amplitude along the x-axis in the four configurations shown in 
Figure 75. The vertical dashed lines show the locations of the point source arrays (magenta) and the inclusion (in 
matching colors). Legend: red - 11x11, green - 12x11, blue - 13x11, orange - 14x11. 

Circular Sonic Crystal 

In light of the above observations on the importance of the space enclosed by the sonic crystal, it is clear 
that a rectangular arrangement has two characteristic lengths, which leads to a higher probability of 
matching the bending wavelength to one of characteristic lengths. For this reason, a different geometry 
is explored: circular arrangements of the sonic crystal. In the simulations, the circular shape is omni- 
directional and has only one characteristic length: the radius (or the diameter) of the ring. 

Two configurations are explored. In both configurations, the outer ring is located at a radius nearest to 
762/2 = 381 mm such that distance between two adjacent inclusions is 57 mm. In the computations, the 
arc length is used to approximate the linear distance. As the result, 42 inclusions are placed at an 
angular distance of 360°/42 = 8.5714° apart at a radius of 381 mm. The first inclusion is placed at the 
polar coordinate (381 mm, 8.5714°/2) such that the upper and lower halves of the space are mirror 
images of each other. The two different configurations come in different arrangements as the inner ring. 
In the first configuration, the inner ring follows the same angular arrangement of the outer ring but is 57 
mm shorter in radius: that is, it has 42 inclusions arranged at a radius of 324 mm. In the second 
configuration, the inclusions are also located at radius of 324 mm, but the distance between adjacent 
inclusions is kept as close to 57 mm as possible. 36 total inclusions are used, placed 10° apart, at an arc 
length of 56.55 mm. Overall, the first configuration maintains the appearance of the "square" grid of 
the sonic crystal; the second configuration maintains the lattice constant as designed. 

The simulation results for both configurations at the target frequency of 3000 Hz are shown in Figure 77. 
Figure 78 shows the corresponding displacement amplitude distribution along the x-axis. In this figure, 
the magenta vertical dashed line represents the location of the point source arrays; and the green 
dashed lines represent the extents of the inclusions in the sonic crystals. 
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The results show that it is still possible to form a configuration resonance. The second configuration, 
which maintains the designed "lattice constant" in both directions, achieves significantly better filtering. 
Although the two configurations have almost identical overall geometrical shapes for the enclosed 
space, the resonance patterns and strengths are very different. 



Figure 77. Displacement amplitude with the presence of circular sonic crystals, at frequency of 3000 Hz. Left: 
maintaining the appearance "square" grid; Right: maintaining the lattice constant as designed. 



Figure 78. Distribution of relative displacement amplitude along the x-axis for two configurations shown in 

Figure 77. 


Concluding Remarks 

These simulations demonstrate the feasibility of using an array of point sources to simulate a planar 
incident wave. Results also indicate that a sonic crystal design based on a single array can perform 
similarly to designs based on planar incident waves that originate infinitely far from the array. In other 
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words, the designed sonic crystal exhibits the bad gap effect regardless of the type of the incident wave 
form. What is important is the frequency of the designed band gap. 

However, when implementing sonic crystal designs in practical finite structures, some practicality 
limitations might significantly degrade performance. In this series of studies, two effects have been 
encountered. The first is due to limited available space to embed the sonic crystal, where only two rings 
of inclusions may be used. This design exhibits the desired band gap but the gap does not act as a total 
blockade. A fraction of the wave can still penetrate the sonic crystal. The second effect is when the sonic 
crystal must be formed as an enclosure around a finite space, where the transmitted waves can form 
standing wave patterns at or near resonance frequencies. In such cases, the filtering effect diminishes 
further. An effective solution to the problem might be to use multiple rows of inclusions. It would be 
worthwhile in the future to explore the effect of having more rows of inclusions, such as using three 
rings of scatterers. 


3.2.5 Split Panel Concept for Airborne Sound Transmission Reduction 

It is well known in the acoustic sound power transmission loss community, and particularly by glass 
window manufacturers, that 'double glazed' panels are preferable to single panels. The same amount of 
mass is spread between two panels with an air gap between them. Using two panels nearly doubles the 
TL when compared to a single panel with the same mass. Figure 79 shows two sandwich panels, with 
thinner 6.35 mm (0.25 inch) cores, with a 12.7 mm (0.5 inch) air gap between them. Each face sheet is 
treated with a center layer of the VHB 9469 discussed in Section 3.2.3. Table 3 lists the thicknesses and 
surface densities of each layer of the optimized panel. 

The 12.7 mm gap is chosen to ensure that sound transmission degradation associated with the mass- 
spring-mass resonance of a double panel system is well below the first transmission tone frequency at 1 
kHz. The resonance frequency, where each panel acts as a lumped mass connected by the stiffness of 
the air gap, is: 


J_ I ~pcFTd_ 

2 n v m x m 2 / ( m x +m 2 ) ' 


(15) 


where pc 2 is the bulk Modulus of air, d is the gap thickness, and mi and m 2 are the two outer panel area 
densities. In the equation, the numerator represents the gap stiffness per unit area, and the 
denominator represents the effective total mass per unit area. The resonance amplifies the sound 
transmission through the double panel system at and around its resonance frequency. The effects of 
the gap thickness on the mass-spring-mass resonance, and on the overall panel thickness, are 
summarized in Table 4. A 12.7 mm gap shifts the resonance below 500 Hz, which is sufficiently below 
the 1 kHz target so that degradation should not occur. 

Rather than leave the air gap empty, Bell suggested filling it with a 9.5 mm (0.375 inch) thick layer of 
Amber Microlite AA insulation (24 kg/m 3 ) from Johns Manville. The insulation provides both thermal, as 
well as reduced sound transmission through its added mass. Bell often adds an extra layer of Microlite 
contained within a thin plastic covering on the inside surfaces of its rotorcraft roof panels. However, the 
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layers are costly, and must often be removed when servicing the panels. Including the insulation inside 
the panel is preferable. The added acoustic transmission loss benefits are modest, and due mostly to 
the added mass of the material, as shown in Table 5. In the Table, the 'Mass Law' column corresponds 
to estimated transmission loss due to the mass of the material, computed using the well-known infinite 
panel formulas [6]. 

An initial analytically based estimate of the performance benefits of the split panel damped design, 
using the analytic tools validated against the baseline panel measurements, is shown in Figure 80. The 
optimized panel coincidence dip is higher in frequency, since the split panel cores are half the thickness 
of that of the baseline panel. Also, the face sheet rigidities are lower due to the softer center VHB 
layers. However, the coincidence dip was targeted to lie between the two transmission tones at 1 and 3 
kHz. Future designs may also adjust core thickness and material to shift the coincidence dip to benign 
source locations. The double panel concept nearly doubles TL at both 1 and 3 kHz. Adding the Microlite 
may also improve TL, but it remains to be seen how well it works above the coincidence dip. Finally, the 
effects of sound transmission through the edge composite is not included in this estimate, but will be 
addressed in final calculations compared to measurements. 


Face sheet/visco g 
sandwiches 

Honeycomb core 



Honeycomb core y» 

I 

Figure 79: Split panel concept. 
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Table 3. Thicknesses and surface densities of optimized panel 


Material 

Thickness 

Thickness 

Surface 

Surface 


(in) 

(cm) 

density 

density 




(lb/in A 2) 

(kg/m A 2) 

Paint 



0.00006 

0.042 

Face sheet 

0.0079 

0.0201 

0.00044 

0.311 

Vi sco 

0.005 

0.0127 

0.00020 

0.14 

Face sheet 

0.0079 

0.0201 

0.00044 

0.311 

Adhesive 

0 

0.0000 

0.00035 

0.245 

Honeycomb 

0.25 

0.6350 

0.00048 

0.3355 

Adhesive 

0 

0.0000 

0.00035 

0.245 

Face sheet 

0.0079 

0.0201 

0.00044 

0.311 

Visco 

0.005 

0.0127 

0.00020 

0.14 

Face sheet 

0.0079 

0.0201 

0.00044 

0.311 

Adhesive 

0 

0.0000 

0.00035 

0.245 

Microlite (3/8") in 1/2" gap 

0.5 

1.2700 

0.00033 

0.228 

Adhesive 

0 

0.0000 

0.00035 

0.245 

Face sheet 

0.0079 

0.0201 

0.00044 

0.311 

Visco 

0.005 

0.0127 

0.00020 

0.14 

Face sheet 

0.0079 

0.0201 

0.00044 

0.311 

Adhesive 

0 

0.0000 

0.00035 

0.245 

Honeycomb 

0.25 

0.6350 

0.00048 

0.3355 

Adhesive 

0 

0.0000 

0.00035 

0.245 

Face sheet 

0.0079 

0.0201 

0.00044 

0.311 

Visco 

0.005 

0.0127 

0.00020 

0.14 

Face sheet 

0.0079 

0.0201 

0.00044 

0.311 

Paint 



0.00006 

0.042 






Total 

1.083 

2.75 

0.00785 

5.501 

total face sheets 




2.488 

total honeycomb 




0.671 

Total visco 




0.560 

Total adhesive 




1.470 

Total Microlite 




0.228 

Total allowable 

0.625 

1.59 

0.00814 

5.700 


Table 4. Effects of airgap thickness on mass-spring resonance frequency and overall panel thickness of 

optimized panel. 


Air gap thickness (mm) 

Resonance Frequency (Hz) 

Overall panel thickness 
(mm) 

3.18 

911 

18.0 

6.35 

644 

21.2 

12.7 

456 

27.5 
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Table 5. Measured and mass-law based sound transmission loss improvements due to use of Microlite. 


Frequency (Hz) 

dB, Measured by 
vendor 

dB, Mass Law 

500 

2.4 

2.4 

1000 

4.6 

4.7 

2000 

6.6 

6.8 

4000 

8.8 

10.5 


CD 

-O 


c 

o 

C/) 

CO 


c 

JO 



Frequency (Hz) 


Figure 80: Analytic estimates of transmission loss benefits of split panel concept, center panel only. 

3.3 Final Optimized Panel Design 

The final optimized panel balances acoustic performance with structural integrity constraints, as well as 
meeting weight and space goals. The split panel concept is augmented with damped face sheets which 
include embedded VHB viscoelastic material, and filled with MicroLite blankets. The final surface 
density is 5.5 kg/m 2 , which is within the allowable goal of 5.7 kg/m 2 . Schematics of the baseline and 
optimized panel cross sections are shown in Figure 81. Although the optimized panel is thicker than the 
baseline panel, the excess thickness is shifted to outside the fuselage. Bell has confirmed that the extra 
thickness will not affect the transmission or other electrical, mechanical, or hydraulic elements in the 
roof cavity region. 
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Baseline Panel 


Optimized Panel 



Above roof 



Below roof 




Figure 81: Schematics of baseline (left) and optimized (right) panel cross sections (not to scale). 

3.3.1 Structural Modeling 

The optimized panel was modeled using hexahedron, wedge, and tetrahedron finite elements, as shown 
in Figure 82. Each ply of fabric and each layer of VHB were discretely modeled with a single layer of 
elements through the thickness. Each core was modeled with two elements through the thickness. 
Adhesive plies were not included in the model, as they have negligible effect on the structural 
performance of the panel. The panel is symmetric except for an extra ply of VHB at the midplane, so the 
lower half of the panel was modeled and then mirrored with an offset for the VHB ply to create the 
upper half. The uppermost surface of the panel was used to create a two-dimensional mesh of 
quadrilateral elements. These elements were extruded to create hexahedron elements. This ply was 
then copied and translated downward to create the remaining plies in the upper half of the panel. The 
core was created by applying a hexahedron and wedge mesh to geometry imported from the CAD 
software Catia. Once all layers were created, each layer was modified to match the design by removing 
the edgeband for the pre-cured plies and removing the center section and ramping down the edgeband 
plies. To create the ply ramps, the hexahedron elements were split into wedge elements along the 
length of the sides and tetrahedron elements in the corners of the panel. 

Beams, straps and angle brackets, as shown in the top of Figure 83, were used to represent the support 
structure of the 429 roof. The beams were modeled with hexahedron elements, with a single element 
through the thickness of the webs and flanges. The beams were connected with straps across the top 
flanges and angle brackets across the webs. These were each modeled with a single hexahedron 
element through the thickness as well. 

The brackets, straps, beams, and panel are connected with fasteners. These elements are connected 
together in the FEM using CBUSH spring elements in the fastener locations. The bottom of Figure 83 
shows the CBUSH elements, represent by yellow circles. Nodes from connected elements are made 
coincident and then a CBUSH, oriented in the direction of the fastener axis, is used to attach them 
together. Then a PBUSH property is assigned to represent the fastener stiffness, computed using the 
bolt material and area. It is important to note that since these CBUSH elements are connected to nodes 
of solid elements, they will not run properly in NASTRAN SOL 101 (Static Analysis), but they will work in 
NASTRAN SOL 103 (Real Normal Modes Analysis) and was the most efficient way to configure the model 
for dynamic analysis. For static analysis, the panel was analyzed in NASTRAN SOL 101, supported with 
Single Point Constraints (SPCs) at the location of the fasteners that would attach it to the beams. 
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The elements were divided into six material families as shown in Figure 82. The carbon/epoxy (C/Ep) is 
divided into three regions for nominal, bag side, and tool side properties. Plies that are cured with VHB 
or as a solid laminate have nominal carbon/epoxy properties. Plies cured between the core and the pre- 
cured laminates have tool side properties. Plies cured between the core and the bag have bag side 
properties. The nominal smeared properties used for dynamic modeling are similar to those listed in 
Table 2. 

It was observed during the test correlation of the baseline panel that the dynamic response is sensitive 
to the total panel weight. While they are not typically discretely modeled for static analysis, the weights 
of the adhesive, paint, and fasteners are accounted for in the FEM for dynamic performance. A CONM2 
point mass element is added at each fastener location to include its weight. The adhesive and paint 
weights are included by adjusting the density of the elements adjacent to them. Separate material cards 
were created for each modified density section using the mass densities listed in Table 3. The weight of 
the optimized panel assembly was estimated from nominal volume and density of each part. The 
densities were modified so that the weights of each part in the FEM equaled the calculated weights. A 
summary of the weights is shown in Table 6. 


Table 6. Optimized panel weight summary. 


Material 

Weight (kg) 

Weight (lb) 

Carbon/Epoxy face sheets 

3.08 

6.80 

Adhesive 

1.89 

4.18 

Paint 

1.02 

2.24 

VHB 

0.466 

1.03 

Floneycomb Core 

0.422 

0.931 

Stiffeners + paint 

9.46 

20.9 

Clips & Angles 

0.712 

1.57 

Fasteners 

0.373 

0.823 

Total Weight 

17.4 

38.4 
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□ Garbo n/Epoxy Bag Side 

■ Carbon/Epoxy Tool Side 

□ Carbon/Epoxy Nominal 
□Viscoelastic {VHB) 

■ Core 
■Aluminum 



Figure 82: FE model of optimized panel. 
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Straps 




Figure 83: FE model of roof structure. Top - brackets and straps. Bottom - CBUSH springs highlighted in yellow. 
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3.3.2 Structural Integrity Calculations 

While the focus of this project is on the acoustic performance of the panel, it is still necessary to analyze 
the structural integrity of the panel with the critical design loads for the representative roof panel. Skin 
panel strength, ramp strength, edgeband fiber and bearing strength, panel stability and step load 
response have been analyzed. The skin panel (regions 1-4), ramp (region 5), and edgeband (region 6) are 
shown in Figure 84. The ramp is least critical of the three, since the skin panel and edgeband both have 
positive margins (see sections 3. 3. 2.1 and 3. 3. 2. 2) the ramp strength is considered good by comparison. 
The ramp is critical for buckling in the stability analysis in section 3. 3. 2.4. 

3.3.2.1 Skin Panel Strength Analysis 

The skin panel refers to the skin plies over the core. The skin plies can be divided into four different 
regions, each having a different ply layup as shown in Figure 76. The most critical layup is from region 1 
as this area has the least number of plies and less support than the skins in regions 2 through 4. Design 
loads for the core stiffened panel were critical for the upper skin in compression. Upper skin applied 
ultimate loads were based on 150% of limit flight loads and were used to analyze the critical skin region 
using an elevated temperature wet open hole compression allowable. The strength based margin of 
safety for fibers in the x (fwd/aft) direction is +0.2. 

3.3.2.2 Edgeband Strength Analysis 

The edgeband fiber strength analysis is similar to the skin panel strength analysis, except that the 
loading moment is applied to the edgeband instead of being coupled out by the upper and lower panels. 
The strength based margin of safety for this load condition was +0.73 for the fibers in the forward/aft 
direction utilizing an elevated temperature wet open hole compression allowable. 

3.3.2.3 Edgeband Bearing Strength Analysis 

The roof beams are attached to the panel using 0.414 cm (0.163 inch) diameter titanium pins through 
the edgeband. The peak fastener applied ultimate load in the roof panel was 4827 N (1086 Ibf), also for 
the jump take-off load condition. The strength based margin of safety for this ultimate load was +0.19 
for the laminate in bearing. 

3.3.2.4 Stability Analysis 

With the addition of an air gap, the ramp down of the laminate from the core to the edgeband becomes 
unsupported. This was deemed a buckling critical area, so a stability analysis of the panel was 
performed. The panel was held fixed at the points where the panel is fastened to the beams and then 
critical loads and moments as determined from previous analysis were applied. A NASTRAN SOL 105 
(Static Buckling Analysis) solution was computed. The resulting critical buckling eigenvalues are both 
about 1.95, which are substantially greater than 1.0, demonstrating that the six plies in the ramp will 
provide adequate stability under the worst case conditions. 
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3. 3. 2. 5 Step Load Analysis 

The panel has the potential to be stepped on when installed in the aircraft. While this is unlikely to 
occur, it is still necessary to ensure that the panel can withstand this step load. To simplify the analysis, 
only the center portion of the top sandwich panel was considered, as shown in Figure 86. The nodes of 
the upper and lower face sheets are fixed from translating vertically at their supported edges. The panel 
is also fixed in all three translational directions at two nodes. The critical condition for a step load is a 
2670 N (600 Ibf) unit load applied to the center of the panel. This represents a 90.7 kg (200 lb) person 
with a gravitational load of 2.0 and an ultimate to limit factor of 1.5. The step load is applied as a 
pressure load distributed across a 12.7 cm (5 inch) diameter area on the upper surface of the panel. 

A NASTRAN SOL106 (Non-Linear Static Analysis) solution was run with this condition. The analysis 
results show that the maximum deflection of the panel is 2.72 cm (1.07 inches), as shown in Figure 87. 
This means that the top panel will contact the bottom panel (since the 1.27 cm [0.5 inch] gap is less than 
2.72 cm [1.07 inch] deflection), and they will share the load. This simplified analysis assumes the top 
panel carries the entire load, therefore the results are conservative. 

The stresses at the center point of the panel for the limiting 2.72 cm (1.07 inch) deflection are: 

a x — —5296 psi 
o y = —8374 psi 
t xy = —10.27 psi 

The maximum forces in the upper face sheet were calculated from the local stresses as shown. The face 
sheet thickness is 0.0079 inches. 

N x = a x t = (—5296 psi)(0.0079 in) = -41.84 

N y = a y t = (-8373 psi)(0.0079 in) = -66.15 

N xy = t xy t = (-10.27 psi)(0.0079 in) = -0.0812 ^ 

These loads were used to calculate the critical strain in the upper face sheet, Ecrit = -907 ps. The 
allowable strain for the face sheet materials is -3575 pe. The positive margin of safety calculated below 
indicates the face sheet is of adequate strength for the step load condition. 

MS = - 1 = 2.94. 

E crit 
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Figure 84: Upper skin plies. 
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Figure 85: Fasteners through roof beams and edgeband. 



Figure 86: Panel section for step load analysis. 



Figure 87: Non-linear analysis deflection results. 
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3. 3. 2. 6 Optimized Panel Manufacturing 

The optimized panel was made from two subpanels, joined together along the outer edge with VHB, and 
via fasteners to the beam frame. Figure 88 - Figure 92 summarize the manufacturing process at Bell- 
Textron. First, an Aluminum tool was machined to facilitate the subpanel layups. Layers of pre-cured 
fabric and VFIB were laid on the flat portion of the tool, and vacuum compacted between layups. 
Uncured fabric was laid in the transition region and edges, and hand worked to ensure proper shaping. 
The honeycomb core was attached, and the rest of the plies overlaid. Each subpanel was then vacuum 
sealed and cured in an autoclave. The two panels were assembled with VFIB, which is also an adhesive, 
and the frame was overlaid to match-drill holes for the fasteners. The fasteners were attached and 
checked for tightness. Finally, holes were drilled along the edges to match those in the NASA SALT 
facility. The final optimized panel is shown in Figure 93. 




Figure 88. Aluminum tooling (top), with layers of fabric and VHB applied (bottom). 
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Figure 89. Top - hand working of edge material. Bottom - application of honeycomb core. 
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Figure 90. Top - Completed single subpanel, Bottom - Curing of vacuum-sealed subpanel. 
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Figure 91. Match drilled holes for framing. 
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Figure 92. Top - tightness check of frame fasteners, Bottom - drilling of mounting holes for NASA SALT facility. 
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Figure 93. Final optimized panel - top - interior face, bottom - exterior face (with frame). 
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4 Optimized Panel Performance 


4.1 Measured Modes, Loss Factors, and Mobilities 

Modal and sound power measurements for the optimized panel were made at ARL/Penn State and in 
NASA's SALT facility. A rectangular grid of measurement points with 2 inch spacing was used for the 
modal measurements. Modes were measured for free and mounted conditions. Samples of the 
extracted modes are shown in Figure 94. 

Resonance frequencies and loss factors for the free panel measurements are shown in Figure 95. The 
modes are split into those dominated by the panel edges and the center panel region. Damping in the 
edge modes is high, ranging between 0.15 and 0.30, while damping in the center panel region ranges 
between 0.04 and 0.07, similar to that in the test coupon with treated face sheets. Recall the baseline 
panel damping without embedded VHB is about 0.01. Modal loss factors of the center panel modes 
measured in the free and mounted panel are also compared in Figure 95. The center panel modal 
frequencies do not change significantly when the panel is mounted, but mounting in SALT does increase 
loss factors slightly. The edge modes, of course, are much stiffer when the panel is mounted, and are 
not readily extracted from the data due to the extremely high loss factors. 

In fact, the damping is so high that all modes become difficult to extract and identify above 500 Hz. 

Edge, center, and beam modes are strongly cross-coupled and the panel response due to drives on the 
center panel and edge regions becomes nearly constant with frequency, as shown in Figure 96. The 
nearly constant response indicates high modal overlap, caused by high damping as well as high modal 
density. Figure 97 compares the averaged center panel transverse vibration response, averaged over 
three drive locations. The figure shows that the optimized back panel response is several dB lower than 
that of the driven front panel above 500 Hz, with reductions increasing with increasing frequency. The 
subpanels are therefore well isolated from each other, which should lead to reduced noise transmission. 
The driven center panel response in the optimized panel is comparable to that of the baseline panel 
below 2 kHz, and mitigated above 2 kHz. The optimized panel response is clearly more highly damped 
than that of the baseline panel. 

Note that the response in Figure 96 and Figure 97 rolls off significantly with increasing frequency. This 
may be due to low input forces from the shaker. Mobilities above 3 kHz may therefore be biased low. 
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Figure 94. Mode shapes of optimized panel. Upper left (1,1); upper right (2,1); lower left (1,3), lower right - 

edge mode. 
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Figure 95. Modal loss factors of optimized panel, free and fixed (in SALT) edge boundary conditions. 



Frequency (Hz) 

Figure 96. Surface-averaged transverse v/F for center and edge regions of optimized panel with free boundaries. 
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Frequency (Hz) 



Frequency (Hz) 

Figure 97. Surface averaged center panel transverse v/F averaged over three center panel drives, optimized 
(front and back center panels) and baseline panels with free boundaries. Top - narrow-band. Bottom - one- 

third octave band. 
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4.2 Radiated Sound Power 

Both structural and acoustic (diffuse field) drives were applied to the panel while mounted in the SALT 
facility (Figure 98 and Figure 99). Transverse forces drove the panel at each of the four joints 
(representing the mounting locations of the transmission) to measure drive point mobility, shown in 
Figure 100. The shaker drive data is useful up to 5 kHz. Given the similarity between the mobilities 
(with the exception of joint 2, which was likely measured at a slightly different location along the beam 
flange), radiated sound power was measured for a drive on a single joint (joint 4). Radiated sound 
power for the structural and diffuse field excitations was computed by summing over a square array of 
intensity measurements. The convergence study for a diffuse drive in Figure 101 shows that a 10.16 cm 
(4 inch) spacing, corresponding to a 13 x 13 array of measured intensities, provides accuracy within 0.1 
dB. The convergence study also shows the presence of the panel-cavity-panel resonance near 400 Hz. 
However, since the SALT facility Schroeder frequency is about 350 Hz, subsequent plots are limited to 
frequencies from 315 Hz upward. 



Figure 98. Optimized panel (inner surface shown) mounted in NASA SALT facility. Reference accelerometers are 
placed on center and edge panel regions. Intensity probe array is to the left. 
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Figure 99. Shaker drive (via stinger and impedance head) on upper left frame corner. 

The sound radiated from the baseline and optimized panels for a transverse joint drive is compared in 
Figure 102. The sound radiated by the optimized panel is much lower than that from the baseline panel, 
with improvements ranging from 6 to 20 dB, with the highest improvements above 2 kHz. The reduced 
sound power for structure-borne drives is due to the significant added damping in the panel, as well as 
damping sandwiched between the beams and panel. 

The sound power transmission loss (power input by the diffuse field divided by the power radiated by 
the panel) for the baseline and optimized panels is shown in Figure 103, and the noise reduction for 
both structural and acoustic drives is compared in Figure 104. Once again, the sound power radiated by 
the optimized panel is much lower than that of the baseline panel, with 6-12 dB improvements. The 
frequency with the smallest improvement is the coincidence dip at 2500 Hz, which was positioned by 
design to be between the 1 and 3 kHz gear mesh tones of highest interest. The TL improvements are 
due to a combination of the split panel design, as well as the increased structural damping and MicroLite 
blanket in the air cavity. The improvements are not as high as anticipated prior to the testing, however. 
Analytic modeling of the transmission loss is therefore used to analyze the performance in more detail. 
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Figure 100. Drive point mobilities at each of the four beam joints. 



Figure 101. Convergence of sound power measurements for diffuse acoustic field drive on optimized panel in 
NASA SALT facility with increasing intensity scan resolution. 
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Frequency (Hz) 

Figure 102. Radiated sound power from a drive at a single joint, optimized vs. baseline panel. 



Frequency (Hz) 

Figure 103. Sound power transmission loss, diffuse field excitation, for optimized and baseline panels. 
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Figure 104. Noise reduction for structural and acoustic (diffuse field) drives. 

4.3 Analytic Diffuse Field Panel TL Modeling 

The high panel damping and resulting high modal overlap at frequencies above 500 Hz allows the use of 
simple analytic models to estimate sound power transmission loss. This simplification results in 
significant computational savings, both in model construction and computer analysis. Infinite panel 
transmission loss theory is therefore used for both the center sandwich panel region, as well as the edge 
region to estimate transmitted sound power. 

Incident power is computed simply as the product of the surface incident intensity and the panel region 
areas (30 in x 36 in = 1080 in 2 for the center panel, and 2 x 5 in x 46 in + 2 x 3 in x 36 in = 676 in 2 for the 
edge regions; see Figure 6 for panel dimensions and subtract a one inch wide frame around the panel 
when clamped in the NASA SALT facility). 

The edge region of the optimized panel is constructed of 12 layers of 0.0079 inch thick carbon fiber 
sheets (the same materials used for the center panel), sandwiched around a single layer of the 0.005 
inch thick VHB material. Each set of six layers above and below the VHB is oriented to give nearly quasi- 
isotropic elastic behavior. The total edge thickness is 0.1 inches. The equivalent Young's modulus of the 
optimized panel edge region is assumed to be 80% that of the baseline panel due to the flexibility of the 
VHB. The measured loss factor of 0.30 is applied to the optimized edge panel analytic model, and the 
lower bound of the measured center panel loss factors - 0.04 - is applied to the center panel. 

The transmitted powers are summed, and combined with the total incident power to compute overall 
panel TL. Figure 105 shows the incident power, power transmitted through the center panel region, and 
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power transmitted through the edge panel regions for the optimized panel. No correction is made for 
the mass-spring-mass resonance of the center panel, which occurs near 500 Hz. The added transmission 
loss of the MicroLite is included in the center panel radiated power. Due to the improved transmission 
loss of the center panel region, the power radiated through the edge is dominant over most of the 
frequency range, with the exception of the coincidence frequency of the center panel near 2 kHz. 

Figure 106 compares the simulated TL for the center panel regions to those of the overall panel TL. 
Clearly, the overall panel TL improvements are not as significant as those of the center panel regions 
alone. Figure 107 compares the analytic total sound power transmission loss for the optimized and 
baseline panels to those measured in the SALT. The agreement is reasonable between 400 and 4 kHz - 
within 3 dB. The center panel coincidence dip is sharper in the measured data, suggesting that the 
MicroLite gap fill material may not be providing the anticipated additional TL reductions. The measured 
edge region coincidence dip in the optimized panel seems to actually be near 8 kHz, instead of near 5 
kHz in the analytic model. This difference distorts the analytic TL curve above 5 kHz, worsening the 
agreement between simulation and measurement. 

The figures indicate that the optimized panel TL is limited by transmission through the edge panels, such 
that the large increases in TL in the center panel do not fully translate to comparable improvements in 
the overall panel. This is particularly true of the MicroLite, which has only a small impact on total panel 
TL near coincidence (and which is not observed in the measured data). Future design improvements to 
roof panels should therefore be complemented by accompanying noise control studies on other 
fuselage panels, which may cause flanking noise transmission paths which are stronger than those 
through the roof panel. 

To investigate the impact of the added damping in the center and edge panels, the optimized panel 
analytic TL is recomputed using the lower damping of the baseline panel. The results in Figure 108 show 
that, consistent with panel theory, added damping increases random incidence TL at and above the 
center and edge panel coincidence dips. The improved TL at 1 kHz is due solely to the split panel 
concept, whereas the TL at 3 kHz is slightly enhanced by the added damping. Above 3 kHz, the edge 
panel damping benefits are more significant. 

In the future, higher fidelity simulations of both TL, as well as sound power induced by structural drives, 
may be made using FE/BE (which is expensive), or SEA modeling (which is more economical). The 
beams, center panel, and edge panel regions may be represented by statistical models of mode types, 
and coupling coefficients to include the effects of the beams in the TL calculations. 
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Figure 105. Incident and transmitted powers for optimized panel. 
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Figure 106. Transmission loss for baseline and optimized panels. 
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Figure 107. TL of baseline and optimized (including MicroLite) panels, simulations and measurements. 
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Figure 108. Analytic TL of baseline and optimized (including MicroLite) panels with variable damping. 
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5 Summary and Conclusions 

Baseline and optimized composite roof panels representative of those used on Bell helicopters have 
been designed, constructed, analyzed, and tested. A structural-acoustic model of the baseline panel has 
been created and analyzed using finite element and boundary element models. Modes of vibration, 
forced vibration response, and sound power transmission loss (TL) all compare well with corresponding 
measurements made on the physical panel in NASA's SALT facility. The numerical modeling approach, 
or Virtual Transmission Loss (VTL) modeling, has therefore been confirmed using the baseline panel 
data. Analytic models of the TL also compare well to the measured and numerical data, showing that 
the simpler analytic approaches may be used to design and assess the optimized panel. 

The analysis and testing results show that the center sandwich panel region of the baseline panel is 
extremely stiff and lightweight, with bending waves supersonic above the coincidence frequency of 
about 700 Hz. Since the transmission tones of interest are between 1 and 4 kHz, this means the baseline 
panel center region would reradiate those tones perfectly. To reduce radiated sound, we have reduced 
structural vibration with added damping, and subdivided the panel into thinner panels separated by an 
airgap to reduce airborne sound transmission. The final optimized panel has face sheets with embedded 
layers of 3M VHB material, increasing structural damping in the center sandwich panel from about 0.01 
to 0.05, and two subpanels, nearly doubling sound power transmission loss through the center sandwich 
panel region. The optimized panel mass density is less than that of the baseline panel with added 
acoustic treatments, but is slightly thicker. To accommodate the extra thickness, the panel is offset 
slightly upward toward the transmission to maintain the same cabin roof height. Several structural 
integrity assessments were conducted on the optimized panel to ensure its suitability for use in a 
commercial rotorcraft. 

Analytic studies of the potential vibration reduction benefits of arrays of inclusions in infinite panels 
were promising, but did not translate well to finite panels with low wavenumber vibrations. The 
inclusion arrays were therefore not manufactured in the optimized panel. The analytic procedures, 
however, have been updated to allow for finite panel effects, and may be used in the future to pursue 
other more effective potential noise reduction designs. 

The optimized panel was tested in the NASA SALT facility, and shows structure-borne noise 
improvements of 6-20 dB, and airborne noise improvements of 6-12 dB. Analytic transmission loss 
studies agree well with the measurements, and show that the edge paneling dominates sound power 
transmission in the optimized panel, limiting the benefits of the center panel quieting measures. 

Several key lessons learned from this study are listed below. 

• Material properties of sandwich panels are not always well understood, even when provided by 
suppliers of the component materials. Constructing small test coupons, measuring structural 
mode shapes and frequencies, and updating material properties using a corresponding FE model 
is therefore essential. 
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• When designing roof panels (or other panels), the thicknesses and material properties should be 
adjusted to tune coincidence dips to lie between key transmission noise tonal frequencies. 

• Highly damped structures with moderately high modal densities may be assessed using simple 
infinite panel transmission loss theory, at significant computational and manpower savings 
compared to FE/BE modeling. 

• Minimizing the sound radiated by one component of a complex multi-component structure does 
not translate to similar total structure noise reduction. In the case of the roof panel investigated 
here, the optimized panel TL is dominated by edge panel sound radiation. Including additional 
MicroLite, therefore, had little added noise reduction benefit. In future designs, all components 
must be considered simultaneously when minimizing sound radiation. 
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Appendix A - Baseline and Optimized Panel Testing 

Structural acoustics measurements of the baseline and optimized panels were made to help validate the 
numerical models, and confirm improvements in structure-borne and air-borne sound power 
transmission loss. The free-free modal analyses were conducted at Penn State ARL and transmission 
loss testing was performed in the SALT facility at NASA Langley. 

A.1 Modal Analysis 

Experimental modal analysis was performed on the panel with free boundary conditions, and installed in 
the TL window. The results from these tests included natural frequencies, mode shapes, loss factors, 
energy, conductance, radiated sound power and estimates of material stiffness (Ex, Ey). The effects of 
the ribs on the vibration were quantified and the modeling procedure for the panel was verified and 
refined based on the test results. The free boundary condition tests were used to confirm the panel 
modeling, and the installed tests used to confirm the boundary condition simulations (ideally clamped, 
but more realistically slightly different from purely clamped). 

The modal test was performed using the roaming hammer approach with six reference accelerometers. 
One of the accelerometer locations was at a vertical location consistent with previous Bell flight 
measurements of a Bell 429. A grid of drive points with approximate spacing of 2" was laid up on the 
flat side of the panel. The grid points were refined near the I-beams and cross-beams so that there are 3 
points across the beam. This required approximately 600 drive points on the panel surface. Additional 
drive points were specified on the beam surfaces, after examination of computational modal analysis 
results from the FE model of the baseline panel. The accelerometers were placed at various heights and 
widths along the panel with at least one on/near each I-beam. Three averages were obtained at each 
drive point to reduce measurement noise. Point mobility and conductance measurements were 
obtained at the reference accelerometer locations, as well as at each of the four beam intersection 
points using a refined procedure with 16 averages which circle around the accelerometer. For the 
modal test conducted while installed in the TL suite, sound intensity sweeps were also performed for 
shaker drives at the reference accelerometer locations (see below). 

The hammer and accelerometer data were collected using PSU in-house software at a sampling 
frequency of 20 kHz for approximately 1.5 sec to achieve a frequency resolution of less than 1 Hz. 
Frequency response and coherence data was checked to ensure data quality. The data were processed 
using PSU in-house modal software which uses singular value decomposition together with rational 
fraction polynomial fitting to produce clean modal data. Time data were also processed using a decay 
method to confirm the modal damping values. 

The free boundary condition test was performed at Penn State prior to arrival at the SALT facility. At 
SALT, the set up time was 4-8 hrs and the estimated measurement time was 8-16 hrs, including 
background noise measurements. The TL tests were performed twice, once before the Transmission 
Loss tests, and again afterward (removing and re-installing all sensors) to establish repeatability 
uncertainty. Several hammer excitations were also repeated. 
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List of hardware (PSU) 


• Force hammer (PCB 086CO2, 50 mV/lbf) 

• Six accelerometers (PCB 352-C67, 100 mV/g) 

• Data acquisition (Nl PXI w/ PXI 8106 controller & 2 PXI 4496 16 channel cards) 

• Monitor, keyboard & mouse 

• Glue and accelerator 

• BNC cables and connectors 

• Chalk line/marker 

List required hardware (SALT) 

• 12 Microphones 

• 5 Intensity probes 

• DAQ 

List required software (PSU) 

• Data acquisition software (PSU developed) 

• Modal analysis software (PSU developed [15]) 

• Matlab 

A.2 Transmission Loss Testing 

Transmission loss measurements were performed on the panels using the intensity technique [16]. The 
panels were mounted with the stiffening beams facing the source (reverberant) room to be prototypical 
of rotorcraft conditions. The reverberant room was excited with an array of loudspeakers and the 
transmitted power were measured using a set of intensity probe measurements on the anechoic side. 
Accelerometers were installed at the same positions as used in the modal test in order to estimate the 
forcing function [17] and correlate TL measurements to the modal tests. The accelerometer data were 
acquired for 1 minute every 10 minutes for ensure adequate averaging and repeatability. 

The frequency range was from 70 Hz - 8 kHz (for the 80 Hz - 12.5 kHz useable octave band frequencies) 
and the sampling frequency was at least 20 kHz. Sound intensity measurements were obtained over a 5' 
x 5' area with a 3" spacing between measurement locations. The tests were conducted with two 
intensity probe spacers to increase the frequency range. Multiple tests were conducted to establish 
repeatability bounds, and with different source levels to establish linearity. The estimated set up time 
was 1-2 hrs and the estimated measurement time was 8 hrs. 

List of hardware (PSU) 

• Six accelerometers (PCB 352-C67, 100 mV/g) 

• Data acquisition (Nl PXI w/ PXI 8106 controller & 2 PXI 4496 16 channel cards) 

• Monitor, keyboard & mouse 
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• Glue and accelerator 

• BNC cables and connectors 

List of hardware (SALT) 

• Five intensity probe kits (B&K 3595, consisting of a Vi" microphone pair (B&K 4197) dual 
preamplifer (B&K 2683) and either a 8.5 mm, 12 mm, or 50 mm spacers) 

• Intensity probe traverse controlled by LabVIEW 

• Sound intensity calibrator (B&K 4297) 

• Three sound intensity conditioning amplifiers (B&K Nexus 2691 4 channel conditioning 
amplifiers) 

• Twelve microphones (GRAS Type 40AQ. > 2 -inch prepolarized, random-incidence microphones 
with Type 26CA preamplifiers) 

• Signal conditioner (PCB 584A 16-channel signal conditioner) 

• Switch board (Precision Filters Inc. 464k switching system) 

• Data acquisition (Nl PXI-1045 chassis with 5 PXI 4472B 24-Bit 8 channel cards) 

• Eighteen high frequency sound sources (Three JBL 2446H and Fifteen BMS 2" compression 
drivers) 

• Six low-frequency sound sources (JBL JRX 1151 15-inch two-way speakers) 

• Four multichannel amplifiers (RANE MA6S amplifiers) 

List of software ( SALT) 

• MATLAB Data Acquistion Toolbox 

• Custom MATLAB GUI 

• Microsoft Excel 

• Nl LabVIEW 

A. 3 Shaker Testing 

A vibration shaker was used to excite the structure. First, a random noise excitation was used with the 
shaker mounted on the bracket. The shaker was hung perpendicular to the panel from a PSU fixture. 
Intensity sweeps were performed to determine the sound radiation. 

List of hardware (PSU) 

• Six accelerometers (PCB 352-C67, 100 mV/g) 

• Data acquisition (Nl PXI w/ PXI 8106 controller & 2 PXI 4496 16 channel cards) 

• Monitor, keyboard & mouse 

• Glue and accelerator 

• BNC cables and connectors 

• Shaker (Wilcoxon F4, 100 mV/g, 100 mV/lb) 

• Shaker mounting fixture 
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Appendix B - Virtual Transmission Loss for an Aluminum Panel 


The sound power transmission loss of an aluminum panel was simulated using a modal approach. The 
physical and material properties used for the panel are listed in Table B-l. The panel was modeled using 
finite elements with simply-supported boundary conditions. 1225 linear plate elements were used. The 
normal modes were computed up to 3000 Hz using NX NASTRAN (93 total modes). The FE mesh is 
shown with a few representative modes in Fig. B-l. A one-to-one boundary element mesh of the plate 
was also created in the baffled condition. The critical frequency of the panel is 1212 Hz. 

Table B-l. Properties of the aluminum panel 


Dimensions 

1.0 x 1.0x0.01 m 

Loss Factor 

0.01 

Density 

2700 kg/m A 3 

Modulus of Elasticity 

69e9 GPa 


To simulate a transmission loss measurement, the panel was excited with a virtual diffuse field pressure. 
Since the spatial correlation of a perfectly diffuse field is a sine function [11], the pressure cross-spectral 
density matrix of the forcing function can be written as 


G FF - 


in ( k o4j) 


sin 

G » G " My 


(B-l) 


where A ,y is the separation distance between points i and j, k 0 is the acoustic wavenumber and G pi is the 
power spectral density (PSD) at point i. Assuming unit PSDs at all points, the sound radiation due to the 
diffuse field was calculated using the CHAMP approach [3]. 

To compute a transmission coefficient, the power incident on the panel was estimated. For a perfectly 
diffuse field in a room, the power incident on an area S is defined as 


CW- 

P inc=^ L S, (B-2) 

where w,„ is the reverberant energy density and c is the sound speed [3]. Using the blocked pressure 
assumption we can approximate w in as: 


W- 

in 



(B-3) 


where p is the acoustic pressure at the boundary and p 0 is the fluid density. Since we apply a unit 
pressure loading to the panel, the squared pressure in Eq. B-3 must be unity. The transmission 
coefficient then becomes 
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T = 


(B-4) 


p 

rad 

S/Ap Q c 


and the virtual transmission loss (VTL) becomes 


VTL = 10 log 


f S/4p 0 c' 

P 

V rad J 


(B-5) 
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Figure B-l. The finite element mesh of the simply supported aluminum panel with 3 representative 

modes. 

As a theoretical comparison, we computed the oblique angle of incidence transmission coefficient for an 
infinite panel defined as 


,co) = 


[2 A c/sin ^] 2 


[ 2 A c / sin0+(D/a>)ri(k o sin^) 4 J +^a>ph-(D/ co)(k 0 sin 


(B-6) 


where co is the angular frequency, cp is the angle of incidence, D is the structural rigidity, rj is the 
structural loss factor and ph is the plate surface density [6]. For honeycomb sandwich panels, an 
effective thin panel D may be inferred from the effective wavespeed described in 2.6.1. The diffuse field 
transmission coefficient is then found using 
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(B-7) 


TZ/I 

| r(^,<w)sin^cos^^ 

M = 

| sin (/> cos (f>d(f) 

o 

Figure B-2 shows the VTL compared to the random incident transmission loss for an infinite panel as well 
as the theoretical mass-law up to 2 kHz. The virtual TL curve captures the coincidence dip well and 
generally follows the mass-control law at low frequencies and the stiffness-control law at high 
frequencies. Due to the low modal overlap at low frequencies, dips in the VTL are seen at the natural 
frequencies. Fig B-3 shows the one-third octave band VTL. To convert VTL to one-third octave band 
frequencies, the mean transmission coefficient within each band was calculated. This is equivalent to 
computing the OTO band spectrum for the radiated power spectrum and the incident power spectrum 
and taking the ratio. 


71 jl 

= | r ($),&>) sin 

o 



Frequency (Hz) 


Figure B-2. Virtual transmission loss for the aluminum panel shown with the theoretical mass law and 
random-incidence transmission loss curve for an infinite panel with the same material properties. The 

value at the critical frequency is shown as a red circle. 
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Appendix C - Determination of Panel Wavespeed and Young's Modulus 


C.l Theory 

The Young's modulus ( E) of a material partially determines the bending wavespeed (cb) of that 
structure. Using the expression for wavespeed of a thin or thick plate, E can then be computed 
once Cb is known or estimated. The following section describes an approach to estimate Cb based 
on the wavenumber transform. 

Shepherd and Hambric recently showed the modal wavenumber (k m ) of a thin, rectangular plate 
with simply supported boundaries can be found using wavenumber transform [12]. After taking 
the one-dimensional wavenumber transform of the plate's mode-shape in the x or y direction, the 
modal wavenumber is found using 

k m =k zl -2n/L, (C.l) 

where k z i is the wavenumber of the first zero after the peak in the wavenumber spectrum and L is 
the length of the plate in the dimension of interest. The modal wavenumber can equivalently be 
determined by adding 2k I L to the wavenumber at the first zero preceding the peak. This 
relationship is shown graphically in Fig. C.l for mode orders 1 and 2. As shown in Eq. C.l, the 
accuracy of the wavespeed estimate then depends on the ability to compute accurately the zeros 
in the wavenumber spectrum. 




Fig. C.l. The wavenumber spectrum for the first- and second-order modes of a simply-supported 
rectangular plate. The spacing between zeros 2k I L can be used to determine the modal wavenumber by 
finding the first zero after the peak wavenumber and subtracting 2k / L . 


Since the modal wavenumbers in each direction are estimated separately, the total wavenumber 
is taken as the square root of the sum of the squares, k = Jkl + k 2 y . The modal wavenumber is 
associated with the natural frequency of the mode so that the wavespeed at that frequency is 
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easily computed as c b = (o n / k n , where co n is the n th natural frequency. The Young's modulus 
can then be estimated using thin plate theory: 


£ = 12p(1 p^JL , (C.2) 

where p represents the material density and h is the plate thickness. If the plate thickness is large 
enough such that shear effects and rotary inertia cannot be neglected, thick plate theory should be 
used. For non-isotropic materials, an estimate of E using Eq. C.2 would represent an effective 
modulus. The independent estimates of k x and k y could also be used to estimate effective 
orthotropic material properties. 


C.2 Implementation using measured mode shapes 

For most practical problems, structural mode shapes are measured using spatially discrete data 
such that the discrete Fourier transform (DFT) must be used for wavenumber analysis. Using the 
analogy from time domain signal processing, the wavenumber resolution will be defined by the 
length of the plate, dk = 2n / L . Since many practical structures are smaller than 2n « 6.28 , the 
wavenumber resolution will be poor without introducing zero padding. The wavenumber 
resolution will thus be increased by the zeropad factor z p [13]. Additionally, the maximum 
wavenumber will be set by the spatial sampling and is related to spatial aliasing of the 
modeshape. An equivalent Nyquist wavenumber can be defined as k nyq = n/dx , where dx is the 

spacing between known vibration points in a coordinate direction. Therefore, for any given 
structure, the wavenumber resolution can be arbitrarily reduced while the maximum 
wavenumber is fixed. 

Mode shapes are best estimated using experimental modal analysis techniques, which typically 
involves acquiring acceleration to force transfer functions at points distributed of the structure 
and computing the complex mode indicator function (CMIF) [14]. Optical methods could also 
be used to estimate the mode shapes after which the wavenumber transform is used to estimate 
the wavespeed for each mode. The individual wavespeed estimates can then be used together to 
perform a least-squares fit of the data such that E is estimated by minimizing the error over all 
the wavenumber estimates. Simply supported boundary conditions are difficult to create in 
practice so that exact modal wavenumber estimates cannot be achieved. However, mode shapes 
for plates with free or clamped boundaries anecdotally to follow similar trends such that the 
reasonably accurate wavespeed estimates are achieved using Eq. C. 1 as long as the mode shape 
data is accurate. Additionally, using a least-squares fit procedure can reduce the biases of 
individual modal wavenumber estimates. 

Wavespeeds were estimated using the discrete wavenumber transform for a 1.0 x 0.7 x 0.005 m 
simply-supported aluminum panel (£=69 GPa, p=2700 kg/m 3 ) up to 1 kHz. Analytical mode 
shapes up to 1 kHz were discretized with 41 points in the x direction and 29 points in they 
direction. The Nyquist wavenumber was 125.7 1/m with a zeropad factor of 20, corresponding to 
a wavenumber resolution of 0.3142 1/m. The structural wavelength at 1 kHz is 0.2204 m which 
equates to nearly 9 points per wavelength (ppw). 
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Figure C.2 shows the exact wavespeed (using thin plate theory) with the wavespeeds inferred 
using the wavenumber transforms. To determine the effect of the discretization on wavespeed 
estimates, the modeshapes were downsampled by a factor of two (dx =. 05 m, ppw=4.4) and four 
(dx= 0. 1 m, ppw=2.2) and the wavespeed estimates were repeated. All estimates show good 
accuracy indicating that the Nyquist criteria of 2 points per wavelength is acceptable for inferring 
wavespeeds using wavenumber transforms. Additionally, the estimates were repeated with a 
reduced zeropad factor of five (dk= 1.26) as shown in Figure C.3. The error in wavespeed 
estimate increases since the resolution in wavenumber decreases and the nulls used to infer the 
wavenumbers are not well resolved. This suggests that minimum zeropad factor is required to 
achieve the appropriate level of wavenumber resolution. A zeropad factor is therefore suggested 
of at least 4n/L (dk< 0.5) as a rule of thumb. 



Fig. C.2. Wavespeed estimates of an aluminum plate for different mesh discretizations. Wavenumber 
analysis of the lowest mesh density (2.2 ppw) predicts the wave speed correctly, since the modes are not 
aliased in space. The number of points per wavelength is referenced to the flexural wavelength at 1 kHz. 
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Fig. C.3. Wavespeed estimates of an aluminum plate for wavenumber resolution of 0.34 (z P = 20) and 1.26 
(Z/>=5). The larger dk produces errors in the estimate due to an inability to accurately resolve the zeros of the 
wavenumber spectrum. 
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Appendix D: Evaluation of Diffusivity and Angular Weighting Factors 
using Intensity Measurements 

D.l Introduction 

Theoretical transmission loss (TL) calculations assume that the acoustic field exciting the structure is 
diffuse [6]. The power flow in a diffuse field has equal probability to be incident from any angle, making 
the time-averaged intensity equal to zero. Additionally, the energy density in a diffuse field is spatially 
uniform [18]. These assumptions are nearly true in large rooms with hard walls and low damping at 
frequencies above the Schroeder frequency. 

However, the sound field in these types of rooms is influenced by specular reflections near the walls so 
that the field can no longer be considered diffuse. Waterhouse determined that the actual mean-square 
pressure can significantly vary from the diffuse field mean-square pressure when the distance from the 
boundary is less than A/5 [19]. For analysis at 100Hz and 500Hz, the distance from the wall required to 
obtain diffuse field conditions should be at least 0.68m (26.5 in) and 0.14m (5.5 in) respectively. 

Given these conditions, one may question the accuracy of applying a diffuse sound field excitation to a 
structure for predicting TL and whether an appropriate correction factor can be determined 
experimentally. Another issue is how the presence of an elastic panel or an open window impacts the 
diffusivity. This section details a set of intensity measurements that were made in an effort to examine 
these questions. 

D.2 Intensity Measurements 

Using a 2-microphone intensity probe, the diffusivity of the sound field in the reverberant chamber of 
the NASA LaRC structural-acoustic loads and transmission (SALT) facility was investigated. 

Measurements were made both with the optimized panel installed and with an open aperture. A GRAS 
50AI-C intensity probe was used with a 12 mm spacer. Pressure was recorded for 120 s with the 
standard TL sources turned on (see Fig D.l). 

Measurements were taken from 0° to 180° in approximately 22.5° increments. The 0° locations was 
oriented parallel with the installed panel. Since the probe is located on a thin rod extended out from a 
mounting handle, it was not possible to rotate about the acoustic center of the probe. Instead the 
probe was rotated about its handle. Intensity was acquired at four measurement positions when the 
panel was installed. These locations are listed in Fig. D.2. The measurement location was as close to the 
panel as physically possible but was nonetheless several inches from the panel. 
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Figure D.l. A diagram of structural-acoustic loads and transmission facility (SALT), not drawn to scale. 
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Intensity was computed using both the cross-spectral density (CSD) method (Eq. D.l) and by estimating 
the particle velocity using Euler's equation, the latter method allowing for investigation of both reactive 
and active intensity (Eq. D.2). The active intensity computed using each method was found to be in 
nearly perfect agreement. For the CSD computations, a Hanning window was used with 50% overlap 
and 2 15 points per block. Intensity was then filtered into one-third octave (OTO) bands, a common 
practice which has been shown to reduce statistical variance [20, 10]. 


/</>=- 


Im{G„.(/)} 

(2n)f pAx 


(D.l) 


/„(/) = -iRe{PCT} 


(D.2) 


D.3 Limitations of intensity measurements 

The intensity computed from these measurements is subject to several limitations. First, the diffuse 
field is only partially correlated over space. This will cause the pressure at one microphone to be 
uncorrelated at the second microphone above some frequency. The correlation function of the diffuse 
field pressure between two points is a sine function: 

sin 
kAx 

where Ax is the separation distance (fixed here by the 12 mm spacer) and k is the wavenumber [10]. 

The quantity kAx can be thought of as the scaled ratio of the separation distance to the wavelength. For 
frequencies below 2.4 kHz, the correlation between the microphones is greater than 0.95 (Ax=A/12). 
This occurs when kAx = ti /6 (orAx=A/12) making the measurements only good up to the 2 kHz octave 
band. This is shown graphically in Fig. D.3. 

Another limitation of the intensity measurement is caused by the finite measurement time. The time- 
averaged intensity is zero in a diffuse field only when averaged over infinite time. For finite-time 
measurements, the relative amplitude and apparent direction can still be compared for qualitative 
assessment but will exhibit randomness in amplitude and direction. An additional error is caused by the 
influence of non-normal intensity incident on the microphones. The cardioid directivity pattern of the 
microphones helps to limit this error but does not eliminate it. 

A final limitation comes from measuring intensity in only one-dimension. Energy density cannot be 
measured without knowledge of the particle velocity in the three Cartesian directions and the acoustic 
pressure. Directional energy density, on the other hand, is more appropriate for an intensity probe but 
requires knowledge of the associated solid angle. Such measurements were not practically feasible for 
this setup and time table. 
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Figure D.3. The partial correlation of the diffuse field limits the frequency range of intensity using the 

2-microphone technique. 

D.4 Intensity Results - Panel Installed 

The response of three accelerometers was acquired during the transmission loss testing. The 
accelerometers were located on the reverberant chamber side of the panel and show the panel's 
vibration response during TL testing (Fig. D.4a). The transmission loss is shown in Fig. D.4b for 
reference. 

The intensity diagrams for the OTO bands from 400 Hz to 4 kHz are shown in Figs. D.5-15. The direction 
of the vector indicates the orientation of the intensity probe and the length designates the level 
reference 1 pW/m 2 . The dB scale shows the level for approximate vector length. The measured 
intensity varies over space, angle and frequency partially due to the random nature of the sound field. 



Figure D.4a. The accelerometer response recorded during transmission loss testing, b) The measured 

transmission loss of the optimized panel. 
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Several interesting observations will be given rather than detailing the features of each plot. First, the 
relative amplitudes of the intensity in the middle of the room are small with respect to the intensities 
near the panel for all measurements below 2 kHz. At 2kHz and above, the dominant direction and 
amplitude of the intensity appears to be random at all locations. This is likely caused by the partial 
correlation of the sound field between the microphones. The intensity at the panel shows a dominant 
intensity direction which varies with frequency, sometime switching directions entirely from one OTO 
band to the next. Also, the intensity at the wall (20" from the panel) shows strong power flow at some 
frequencies (e.g. 500 and 800 Hz) but very little at other frequencies (e.g. 400 and 1250 Hz). 
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Figure D.6. Intensity field at 500 Hz with the optimized panel installed. 
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D.5 Intensity Results - Open Aperture 

The panel was then removed from the window and intensity measurements were repeated. 
Measurements at the open window (where the panel would have been located) and ~15" inside the 
anechoic chamber were acquired at all the angles used from the previous measurements. Intensity 
measurements were made only at 0° and 90° for the locations in the reverberation chamber and at 0°, 
45° and 90° for an additional position located approximately 5' inside the anechoic chamber. The results 
are shown for the 1 kHz OTO band in Fig. D.17. As expected, the intensity in the middle of the 
reverberant room is small but acoustic power clearly flows from the reverberant room to the anechoic 
room. A similar trend is seen for all OTO frequency bands. 



Figure D.17. Intensity field at 1 kHz with an open aperture. 

The open aperture measurements were used to estimate an angle-dependent correction factor using a 
procedure similar to that in references [21] and [22]. The normalized intensity was averaged over all the 
OTO bands and is shown in Fig. D.18. The correction factors proposed by references [21] and [22] are: 

cos 12 (6) (D.4) 

g-0.50 

and are also shown in the Figure. 
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The exponential coefficient was modified to 0.5 to provide a better fit than the value of 1.5 used in 
reference [6]. Both correction factors show good agreement from angles from 0° to 45° but deviate 
slightly at larger angles. 



Figure D.17. Normalized intensity averaged over OTO bands with the open aperture for two positions: 
'Panel', which is the open window (where the panel would be) and 'Anechoic', which is inside the 
anechoic chamber. Two possible weighting functions are also shown. 

Simulations for the transmission through a thin aluminum panel (0.813 mm) were then updated using 
the cosine weighting coefficient as shown in Fig. D.18, and compared to measurements made in the 
NASA SALT facility. The simulations were made using the fully integrated diffuse field sound power 
transmission coefficient (DAF theory) and its accompanying low frequency approximate form (DAF 
approx.) as well as the normal incident mass law (equations 5.50, 5.51 and 5.20, respectively from 
reference [6]). The TL estimate is improved above 400 Hz when including the angular weighting. This 
case only illustrates agreement for transmission below the panel's coincidence frequency. 
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Figure D.18. Simulated and measured TLfor an aluminum panel. The simulation is improved when 

using the cosine weighting function. 

It should be noted that using the open aperture intensity field to determine an angular weighting 
function for TL testing has been questioned by Fahy [23], He states that the intensity in a hemi-diffuse 
field incident on a non-reflecting plate inherently has a cos(0) dependence and therefore assigning a 
weighting function for TL measurements from an open window measurement is not physically 
justifiable. This idea is consistent with Pierce's derivation of diffuse field intensity through an open 
window [10]. However, the predicted [21] and measured [22] weighting factors (using beamforming) 
are obtained without the use of an open aperture and are still valid. 

D.6 Future work 

Additional experimental and numerical work could be performed to characterize the diffusivity of the 
reverberation chamber and improve TL simulations. Improved particle velocity measurements should 
be made using 3-D velocity probes to accurately measure generalized energy density at a point, which 
has been shown to have lower spatial variance than pressure measurements [24, 25]. However, the 
expense of such probes is somewhat prohibitive such that other measurement techniques should be 
developed. 

One such procedure could measure pressure at points near a hard wall to map out the spatial 
correlation and statistics of the sound field. The measurement points would vary over the three axes (x, 
y and z) to show a complete representation of the sound field and then be compared to existing theory 
[19,26, 27], The measurements could then be repeated near an elastic structure operating below 
coincidence and compared to those for a hard wall. The measurements would also be used to evaluate 
two proposed weighing functions [21, 22]. Additionally, more advanced metrics for diffusivity could also 
be computed [18, 20]. 


127 


Additional theoretical modeling could also be performed to determine the effect of an elastic panel or 
open window on the sound field in a room. An analytical model could be developed by modifying a 
hybrid Green's function / modal summation method for computing the sound field in a reverberant 
room [25]. The hybrid approach has been shown to improve convergence for rigid-walled rooms. A 
generalized impedance boundary would be included in the model to include the effect of an elastic 
panel or open aperture. Using an analytical model, the sound field could then computed in a 
reverberant room near and away from rigid boundaries to assess the diffusivity and the general impact 
of the incident intensity on TL could be predicted. 


128 



Appendix E - Formulation for Multiple Scattering of Flexural Waves in a 
Thin Plate 

From Dr. Liang-Wu Cai, Kansas State University 
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Formulation of Multiple Scattering Problems 
of Flexural Waves in Thin Plates 


Liang- Wu Cai — Kansas State University 
for the Acoustically Tailored Composite Rotorcraft Fuselage Panels Project 

Last Updated: November 2014 


I. Introduction 

There are several types of waves that can exist in a plate. In this report, only the fl xural 
waves are considered. 

The justificatio for limiting the consideration to fl xural wave is that the primary concern 
for the present project is at lower frequencies, for which the wave length is expected to be 
longer than the plate thickness. For example, the frequency of the primary concern is approx- 
imately 3 kHz. Assuming a wave speed in the order of 1000 m/s, the wavelength would be 
in the order of 0.5 m. Other types of waves, most prominently, the Lamb waves, will become 
important when the wavelength is in the same order as the plate thickness, which is in the 
order of 10 -3 ~ 10 _2 m. 

A good design feature in the plate geometry is that the fl xural wave speed is highly des- 
ignable, although frequency-dependent, by varying the thickness of the plate and in turn its 
rigidity, as it shall become clear in the expression for the wave number. It would be highly 
desirable to reduce the fl xural wave speed so to bring the wave length down to the order of 
0.1 m. Even in this wave length, it is still at least one order of magnitude longer than when the 
Lamb waves become a concern. 


II. Governing Equation and General Solution 

II. 1 General Solution 


The governing equation for a wave propagating in a thin plate can be written as (Pao and 
Mow, 1971, p. 326) 


DV 4 w + ph 


3 2 w 


= q 


(i) 


where w(r,t ) is the deflectio (out-of-plane displacement) of the plate, h and p are the thick- 
ness and the mass density, respectively, of the plate, q is the distributed lateral pressure, and 
D is the plate’s rigidity or flexural stiffness. For classic plate theory, the fl xural stiffness for 
a uniform plate is given as 


D = 


Eh 3 

12(1 - v 2 ) 


( 2 ) 
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where E is the Young’s modulus, and v is Poisson’s ratio. For a composite sandwich panel 
with two identical isotropic face sheets and an isotropic core (Vinson, Plate and Panel Struc- 
tures of Isotropic, Composite and Piezoelectric Materials, Including Sandwich Construction, 
Springer, 2005, p. 296.) 


Efh 2 c tf r \E c i}-vj)h c 

2(1 -vj) |_ 6E f (l -vf)t f 


(3) 


where parameters with a subscript / belong to the face sheets and parameters with subscript 
c belong to the core, tf is the thickness of one of the face sheets. 

With the absence of a lateral loading, the equation becomes homogeneous. In the steady 
state, assuming 


w(r, t) = W(r)e~ ta,t 


(4) 


eqn. (1) becomes 

V*W-—co 2 W = 0 (5) 

which can be split into two equations 

S/ 2 W 1 +k 2 Wi= 0 and V 2 W 2 - k 2 W 2 = 0 (6) 

where k 2 = oj ^ phj D. They are the Helmholtz equation and the Modified Helmholtz equa- 
tion, respectively. In a polar coordinate system, their solutions are linear combinations of 
cylindrical Bessel functions with kr and ikr as respective arguments. Alternatively, the 
Bessel functions with imaginary arguments can be replaced by the modifie Bessel functions. 
Thus, the most general solution in the polar coordinate system can be written as 

OO 

W(r,6) = Y, [AnJnifir) + B n H n (kr ) + C n I n (kr) + D n K n (kr)] <f n6 (7) 

n——oo 

where //„(•) is the Hankel function of the firs kind, /„(•) and K n (f are the modifie Bessel 
function of the firs and the second kinds, respectively. 

Matrix notation is introduced such that a set of wave expansion bases {3f(r, 0)} is define 
as 

i2?(r,Q)} n = {& n (kr)e ne } (8) 

with the index n running from — oo to oo, where S n (') represents any of the Bessel functions. 
Then, the general solution can be written as 

W(r, 6) = {A} T {J(r, 6)} + {B } T {H(r, 6)} + {C} T {I(r, 6)} + {D} T {K(r, 0)} (9) 
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II.2 Modified Bessel Function 


The modifie Bessel functions are define as the following 

In(z) = r n j n (iz ) (10) 

K n (z) = ^t” +1 H™(tz) = ^(-i) n+1 Hf\-iz) (11) 

The general behaviors of these functions are similar to exponential functions: as z increases, 
I n (z ) increases without bound, while K n (z) decreases. Thus, for a problem define in an 
infinit domain, I n (z) would not appear, and K n (z) represents motions confine in the vicinity 
of the origin (z = 0). On the other hand, for a problem whose domain includes the origin, 
K n {z) would not appear, as it is singular at the origin. 

II.3 Moments and Shear Forces Across Thickness of Plate 


The moments and shear forces across the thickness of the plate can be written in terms of 
displacement w as (Pao & Mow, 1971, p 329) 


M r 


M e 


M r g 


Qr 

Qe 

V r 


-D 


vV 2 W + (1 - v) 
D I V 2 W — (1 -v) 


d 2 W' 

dr 2 

d 2 W~\ 


dr 2 


a ndw 


D( l-v)— - 


dr \ r 3 6 


a 




-D 


a 


(V 2 kf) 


rdO 
1 9 M r0 

~r 9(9 


( 12 ) 


where, Q's are the shear forces on the cross-section of the panel, and V r is the so-called total 
(net) shear force that combines the contributions from the shear force and the in-plane moment 
M r g (Graff, Wave Motions in Elastic Solids, Dover, 1975, pp. 234-235). 

Among these quantities, the most important ones are those used in the boundary condi- 
tions. They are M r , Mg and V r . Correspondingly, a set of special functions & s (n,z) is 
defined and the corresponding wave expansion bases {& s {r, 0)} are define as 


{& s (r,9)}„ = & s (n,kr)e ln6 


(13) 


where the superscript t denotes the type of Bessel functions: 1 for J„{z), 2 for H n (z), 3 for 
/„(- ), and 4 for K n (z ); and the subscript ,s' denotes the sequence of moment- force component: 
1 for M r , 2 for Mg, and 3 for V r . Detailed definition for these ©-functions are given in 
Appendix A. 
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( 14 ) 


With this set of ©-function, if a deflectio is expressible as 

W = {ct} T {$r(r,0)} 

The corresponding moments and shear force are expressible as 

M r = -\{ci} T {<£ l (r,6)} 

M e = — \r{a} T {<& 2 (r, 0)} 

r 2 

V r =-^{a} T {<Z 3 (r,9)} (15) 

r 3 

III. Single Scattering Problems 

The single scattering problem lays the foundation for the multiple scattering analysis. The 
analysis of a single scattering problem leads to the establishment of the so-called T -matrix of 
the scatterer; and in the process, the matrix notation is defined 

III.l 7 -Matrix in Generalized Matrix Notation 

Consider a single circular scatterer of radius a that is located at the origin, and subjected 
to a planar incident wave. The incident wave is expressible as 

W' nc = {A P } T {J (r, 9)} + {A e } T {I(r, 0)} (16) 

and the scatterer wave, which exists in the infinit domain, is expressible as 

W sa = {B P } T {H (r, 9)} + {B e } T {K(r, 9)} (17) 

where the superscript p and e in the wave expansion coefficient signify the propagating and 
evanescent components, respectively. 

Physically the {A e } component of the incident wave would not exist in the single scattering 
problem since this component is localized. But in multiple scattering, the wave scattered from 
nearby scatterer would contain this component and acts as a part of the incident wave. 

As with any scattering problems in a linear system, the incident and the scattered waves 
form a cause-result causality relation, and can be represented by a matrix, conventionally 
called the T-matrix. For the problem at hand, there are two components in each of the incident 
and the scattered waves. They form the following set of four relations 

{B pp } = [T PP ]{A P ) {B ep } = [T ep ]{A p ] (18) 

{B pe } = [T pe ]{A e } {B ee } = [T ee ]{A e ) (19) 

where the firs pair of equations determines the scattered wave due to an incident propagating 
wave, and the second pair is due to an incident evanescent wave. In the double superscripts 
for { B } and [T] matrices, the firs superscript signifie the type of the incident wave, and the 
second superscript signifie the type of the scattered wave. 
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In order to maintain the conventional form of the T -matrix relation, the above relations 
can be combined into the following relation in generalized matrices 


{58} = [*]{*} 

where the generalized matrices are define as 


{58} = 


{ \%\ } <«> ’ - { \ A A '\ } m - 


[T pp ] [T pe ] 
[T ep ] [T ee ] 


( 20 ) 


( 21 ) 


Correspondingly, the wave expansion bases can also be concatenated to form the generalized 
wave expansion matrices such that 

{s,r HK!} (22> 

Then, the incident and the scattered waves are representable as 

J^inc = g,)} W scr = 0 )} (23) 

III.2 Rigid Inclusion 

For a rigid inclusion, the boundary conditions are that the plate is clamped at the edge of 
the inclusion. That is W = 0 and dW /dr = 0. This leads to the following equations 

B p H n (ka ) + B e n K n {ka) = -A p J n (ka) - A e n I n (ka ) (24) 

B p H’ n {ka) + B e n K' n (ka) = - A p J'(ka ) - A e n I' n {ka) (25) 

which can be readily solved as a set of linear equations subjected to two different sets of 
right-hand sides. The solutions can be written as 

_ Bj/ P -J n {ka)K' n {ka) + J' n {ka)K n {ka) 


[T pp L 
[T ep ] nn = 
[T pe } nn = 


[T ee ] 


A p 

Sin 

r> e P 

~J T 

sin 

nP e 

“n 

bi_ 


J n {ka)H' n {ka) - J’ n (ka)H n {ka) 21 1 


- I n (ka)K' n (ka ) + I^(ka)K n (ka) 
An 

I n (ka)H' n (ka) - I^(ka)H n (ka) 


nka A n 
1 1 
ka A„ 


where 


A„ = H n {ka)K' n (ka) - H' n {ka) K n {ka) 


and the Wronskian relations of Bessel functions have been used 

Jn(z)K(z ) - J'(z)H n (z ) = — 

71 Z 

I n (z)K' n (z) - r n (z)K n (z) = -- 

Z 


(26) 

(27) 

(28) 

(29) 

(30) 

(31) 

(32) 
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III.3 Void Inclusion 


For a void inclusion, the boundary conditions are that the moment and the force along the 
hole’s edge vanish. That is M r = 0 and V r = 0. This leads to the following equations 


B p <&\(n,ka) + B e n <&\(n,ka ) = — A p &\(n, ka) — A e n <£\(n,ka) 
B p <&\(n,ka) + B e n <$b\(n,ka) = — A^\{n,ka) — A e n <&\(n,ka) 

which can be readily solved to give 

B pp © 3 (n, ka)&\(n, ka) — ©}(n, ka)&^(n, ka) 


[T rP ]„n = 
[T ep )„ = 
\T“\„ = 

^pee 


An 


A^ 


Bn P © 3 (n, ka)&\(n, ka) — ©j (n, ka)&\(n, ka) 

~a[ = a ; 

B pe &l(n, /ra)©j(n, ka) — ©j (n, ka)&*(n, ka) 
A e n A„ 

© 3 (n, ka)&\(n, ka) — ©j (n, ka)&l(n, ka) 
An 


K 


(33) 

(34) 


(35) 

(36) 

(37) 

(38) 


where 


A„ = ©j(n, ka)(£^(n, ka) — ©j(n, ka)&l(n, ka) 


(39) 


III.4 Elastic Inclusion 

For an elastic inclusion, waves exist inside the scattered Note that the scatterer include 
the origin, thus only non-singular Bessel functions can be used in the expression for the wave 
inside the inclusion; that is 

W Kf = {C p } T {J(r, 6)} + {C e } T {I(r,0)} (40) 

The boundary conditions require the continuity of IT, W' , M r and V r . This gives the follow- 
ing set of equations 

A p J n {ka) + A e n I n (ka) + B% H n (ka) + B e n K n (ka) = C p J n (k x a) + C„I n (k x a) 
A p J' n (ka) + A e n I' n (ka) + B p H' n {ka) + B e n K' n {ka) = \ \C P J' n (k x a) + C^I'Jk^a)] 
A%<£;\(n, ka) + A e n &\(ka) + B^\{n,ka) + B e n ^;\{n,ka) = C p &\(n,k x a) + C*&\(n,k\a) 
A p &\(n,ka) + A e n &\(ka) + B p< S:\(n,ka) + B e n $b\{n,ka) = C p< &\(n,k\a) + C^\{n,k x a) 

which can be converted into the following two sets of linear equation systems and solved 
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numerically: 


—H n (ka) -K„(ka) J n (k x a) I n (k x a) 


BP 


J„(ka) 

-H'„(ka) —K' n (ka) J' n (k,a ) ^(fc.a) 

< 

B e n 

> < 

J'Jka) 

—<&^(n,ka) — ©j (n,ka) &\(n,kia) &l(n,kia) 


Cn P 


©}(n, ka ) 

_—&l(n,ka) —<£*(n,ka) © \{n,k\a ) ©3(77, k\a) _ 


c e 


k ©3(77, ka) J 

—H n (ka) —K„(ka) J n kk\a ) I n (k\a ) 


( RP 


' hika ) 1 

—H' n {ka) —K' n {ka) ’fj^a) ’fl^a) 

A 

K 


I' n ika) 

—&^(n,ka) —(£^(n,ka) &\(n,k\a) <&\(n,k\a) 

> 

C P 


©j (n, ka) 

_—&\{n,ka) —©3 (n,ka) ©3(77, Aqa) <&\(n,k\ri) _ 


[cp 


, ©3(77, ka) j 


A' (41) 


A’ (42) 


The solutions will give a set of eight ratios: the ratios between B n and A„ give the entries of 
the T -matrices: 


[T xy \ 


B xy 

Ay 


where v and y is either p ox e. The ratios between C n and A n defin another set of character- 
istic matrices. More often, they are re-expressed with respect to the T -matrices, as 

C xy 


{£} = [ 9 1 ]{ 23 } 


with 


[R Xy \nn - 


B x y 


(43) 


IV. Multiple Scattering Problems 

IV. 1 Generalized Matrix Formulation 

With the introduction of the generalized matrix notation, the general multiple scattering so- 
lution can be readily obtained, following the same solution structure for other types of waves, 
as long as the system is linear. 

In a general problem setup, a set of N scatterers are located within an infinit plate. The 
incident wave can still be expressed as eqn. (16), although in reality, the evanescent component 
would not present. The wave scattered by individual scatterers can be written as 

Wr = { B?} r {H(ri , Oi)} + {Bf} r {K(n, 6, )} (44) 

which remains similar to the one for the single scattering problem, in eqn. (17), with two 
differences: the subscript i is used to signify the parameters belonging to scatterer i, and the 
waves are expressed in local coordinate systems. The “classic” multiple scattering formulation 
gives the following solutions, in terms of generalized wave expansion coefficien matrices: 

{»,} = p-, ]({«,} + Y, [“ufW) (45) 

\ 7 = 1 / 

which can be solved by rearranging into a linear equation system as 

[L]{B} = {A} (46) 
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where 


{B} 


1 



{^ 2 } 

> {&} = < 






(47) 


[L] 


[S'] -[^i][»ti2] r -[‘StI^is] 7 ’ 

-[^ 2 ][^2l] r [»] -[^2][^23] r 




-[^iv][»tiVl] r -[*nW N2 ] T ~[^nWn3] T ••• [»] 


The only remaining issues are matrices [flfcy], which represent the coordinate transformations, 
and matrices {21,-}, which is the wave expansion coefficient for the incident wave. 

IV.2 Coordinate Transformation 

The coordinate transformation among the local coordinate systems is obtained from the 
Graf’s addition theorem for the Bessel functions. It has been a crucial ingredient for the 
multiple scattering formulations, and hence it is well known. However, the modifie Bessel 
functions possess most of essential properties of the Bessel functions and yet differ in some 
properties. Therefore it is useful to trace the derivation process for coordinate transformation 
for the propagating wave expansion bases, and then follow the same process to derive the 
coordinate transformation for the evanescent wave expansion bases. 

IV.2.A Coordinate Transformation for Propagating Wave Expansion Bases 

Recall Graf’s addition theorem for the Bessel functions as (Abramowitz & Stegun, 9.1.79, 
p. 363 1965): 

(X) 

n P = Y Jn+m(Z)J m (z ) °? S nw (49) 

Sill Ijlll 

m=— oo 

and 

(X) 

H n (vT) nf = Y H n+m (Z)J m (z ) °° S ma (50) 

ulll Ijlll 

m = — oc 

where both n and m are integers. Equation (50) is valid only when \Z\ > |ze ±z “|; whereas 
eqn. (49) is valid throughout the plane. When Z, z and m are real and positive, the geometric 
relations of the parameters are as sketched in Fig. 1(a), which is compared with geometric 
relations between two arbitrarily chosen local coordinate systems in Fig. 1(b). An arbitrary 
fiel point P can be located by polar coordinates of either (r t , Of) or (rj , 6 j). Further, (dij , Oij ) 


137 





(b) 


Fig.l Geometries in multiple scattering, (a) Geometry for Graf’s addition theorem for 
Bessel functions. ( b ) Geometric relations between two local coordinate systems. 


is the polar coordinates of Oj in the coordinate system originating at o ; - , and (dji, Oji) is the 
coordinates of Oj in the coordinate system originating at Oj. From Fig. lb 

djj = dji 0(j = TC Oji (51) 

First, eqns. (49) and (50) is rewritten using exponential harmonics as 

OO 

Y n (m)^ = J2 Vn+m(Z)J m (z)e ima (52) 

m=—o o 

where, following Abramowitz & Stegun, Y n (z) denotes either J n (z) or H n (z ). 

Comparing the geometries in Fig. 1, after the following parameter substitutions: 

z — > k r j tu — »• ki'j Z — > kdji a — » Oji — Oj and /3 — 6j — 6 t j, 

eqn. (52) becomes 

OO 

Ynikn)^- 6 ^ = %i+m ( k dj i ) J m (k ij ) & n ( °i ' ~ e i ) (53) 

m = — oo 

Moving Oij from the left-hand side to the right-hand side gives 

OO 

Y n (kri)e ndi = ^ (n0 ’ 1+m0 ’' ) Yn+m(kd ji )J m (kr i )Q- im0 i (54) 

m = — oo 

Since m ranges from — oo to oo, changing — m to m does not change outcome of the summa- 
tion: 

OO 

Y n (kn)e in6 - = J2 eHn6ij ~ mej,) Zn- m (kdji ) J-m (k rj )e im0 J (55) 

m——oo 
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( 56 ) 


Noting that Qp = % — n, and 

Qi(nOij—mOji) __ Qi[(n— m)0ij +mn] __ m )0ij 

and that (A&S, 9.1.5, p. 358) 


(-1 ) m J- m (z) = J m (z) 


eqns. (50) and (49) finall become 


{J(r h e l )} = WR’’ J ]{J(r,,0 J )} (57) 

= (58) 

where the entries of matrices [Rf] and [diRf] at the «-th row and the m-th column are 

[a = e«”-”> e »y„_„(<:rf, 7 ) (59) 

[R!,] nm = H n - m (kd„) (60) 


and the validity condition for eqn. (58) is djj > rj. Matrices [Rfj] and [dtRf] are called the 
singular and regular, respectively, local coordinate translation matrices for the propagating 
wave expansion bases. Here, the symbol denotes the regular counterpart of the enclosed 
function, which is, in this case, obtained by replacing the Hankel function of the firs kind 
by the Bessel function of the firs kind. (Although J n (-) is the real part of H n (-), there is a 
complex factor in the front, it is not appropriate to called [SH Rf] as the real part of [Rfj]- On 
the other hand, H n (•) is singular at the origin, while J n (•) is non-singular.) 

IV.2.B Coordinate Transformation for Evanescent Wave Expansion Bases 

For the modifie Bessel functions, based on the same geometry in Fig. If/), the corre- 
sponding Graf’s addition theorem is (Watson, A Treatise on the Theory of Bessel Functions, 
2nd edn., Cambridge University Press, 1944, §11.3, p3. 61) 

OO 

In{tu) C ° S nfi = (-1 ) m I n+m (Z)I m (z) C ° S ma (61) 

sm sm 

m=—oo 

and 

OO 

K n (m) C ° S nf= J2 K n+m (Z)I m (z ) C ° S ma (62) 

olll ijlll 

m = — oo 

which can be combined to write 

OO 

= Y, (-1 ) m I n+m (Z)I m (z) e ima (63) 

m=—o o 
oo 

K n (vj)e in P = £ K n+m (Z)I m (zW ma (64) 
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Note that a “unified expression similar to eqn. (52) is not available. However, the same 
derivation process can be followed exactly without any change, until reaching eqn. (56). The 
“modified version of eqn. (56) is (A&S, 9.6.6, p. 375) 

I -mi 2 ) = Im(?) (65) 


Thus, 


JmOi 


I n {kri)e indi = ^ e ,(n m)e ^ I n - m (kdij)I m (krj)& 

m=—oo 
oo 

K n (kri)e nd i= J2 (- 1 K n _ m (kd u ) I m (k r-j )e 


imO; 


In matrix form, 


{I(r i ,0 i )} = [MR e ij ]{I(r j ,0 j )} 


{K(r i ,e i )} = [R e ij ]{I(r j ,O j )} 


where 


[diR^ j } nm =^ n - m ^I n - m (kd ij ) 

[Xtfnm = (-\) m e Kn - m)6ii Kn-mikdij) 


( 66 ) 

(67) 


(68) 

(69) 

(70) 

(71) 


IV.2.C Coordinate Transformation in Generalized Matrix Form 

Finally, in generalized matrix form, the coordinate transformation can be written as 


[918t.y] 


0 

0 [SHfiJ] . 



(72) 


Note that there is no off-diagonal blocks. Physically, a propagating wave would not be con- 
verted into an evanescent wave by a mere coordinate transformation, and vise versa. 


V. Incident Wave Fields 

V.l Wave Expansion for Planar Waves 

For the propagating incident wave, a planar wave with unit amplitude propagating in the 
v -direction, is given by the well-known relation 

OO 

e ikx = Q ikr cos 6 = i n J n (kr)e tn 9 (73) 

n——oo 

For the evanescent wave, there is no real meaningful “wave”, rather than a localized vi- 
bration, which will be represented by an exponential decaying function. In particular, t~ k \ x \ 
represents such a vibration localized along the y-axis. 
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According to the following properties of modifie Bessel functions (Abrmowitz & Stegun, 
9.6.38, p.3376) 

OO 

qZ cos o _ j Q ( z j _|_ 2 ^ 7„(z) cos(«$) (74) 

n= 1 


Recall that I- n (z) = /„(z). Further, 

7„(z)e'" 0 + 7_„(z)e - '” 0 = I„(z ) (e'” 0 + e - '” 0 ) = 27„(z) cos(nO) 
Thus, replacing z by kr, 


JtX 


Jkr cos# 


X] In(kr)Q 


jn6 


(75) 


(76) 


Furthermore (Abrmowitz & Stegun, 9.6.30 and setting m = 1, p. 3376), 7„(— z) = (— l)"7„(z) 


e - ** = Q-krcosO = (-l)”/„(A:r)e ?n0 (77) 

n——oo 

The same relation can be obtained from eqn. (16) by replacing k with ik, and utilizing the 
relation between J n (z) and /„(z). 

In terms of wave expansion, if the incident wave is written as 

w mc = {AffiJinA)} + {d e }J {I in, 6i)} (78) 


Then, the coefficient are 


{A P )n = ^ 


{A e }n 


(— 1)” if cos 6 > 0 
1 if cos 9 < 0 


(79) 


Note that in two-dimensional scattering problems in a free space, the incident wave is often 
considered as “located at “infinity” its exact location does not matter while the expression 
remains valid no matter where the origin of the coordinate system is located. However, in a 
plate, this is no longer true. The evanescent component of the wave comes with a location, 
which is the v-axis of the coordinate system. For most wave scattering considerations, this is 
not the main concern and thus be neglected. 

V.2 Coordinate Transformation for Planar Incident Waves 


In multiple scattering formulation, a same incident wave is expressed in different local 
coordinate systems. In general, coordinate transformation utilizing the addition theorems of 
Bessel functions can be used to fin the wave expansion coefficient in different coordinate 
systems. However, it mus be noted in the computation, since the infinit series are truncated 
to finit terms, there are errors introduced into the process. As an alternative, coordinate 
transformation can be directly applied to the analytical expressions before expanding it into 
wave expansions. 
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For example, if the incident wave is expressed in a global coordinate system (without 
subscript), and in the scatterer i ’s local coordinate system, whose origin is located at (A), Yf) 
in the global Cartesian coordinate system and (Rj, 0/ ) the global polar coordinate system, the 
following relation can be observed x = Xj + x* . Then, 

Q?kx __ gik(Xi +Xi) __ gikXj gikxi ^—kx __ Q—kXjQ—ikxi (80) 

Then, the wave expansion coefficient for the propagating component of the incident wave in 
the local coordinate system can be simply written as 

{A?}„ = e kx 'i n (81) 

This is a much more efficien and accurate approach to obtaining the wave expansion coeffi 
cients for the individual scatterers in the multiple scattering problem. However, the expression 
for the evanescent component has not been not sought. 


V.3 Planar Wave Reflected by a Boundary 

At times it is desirable to consider a wave scattering problem in a finit plate. The simplest 
geometry of a finit plate would be a rectangular one. In order to consider the wave reflectio 
from a boundary, the reflecte wave can be approximated as an incident wave propagating in 
the —x direction. As with the conventional incident wave, generally considering the propa- 
gating component is sufficient but the reflecte wave may different by a certain phase. For a 
planar wave of unit amplitude traveling in the — x direction, it can be expressed as e lk( - x °~ x \ 
where x 0 can be a certain reference point and in turn kx 0 = <fi o can be viewed as the phase 
angle. Furthermore, for a reflecte planar, the amplitude of the reflecte wave could be smaller 
than the incident; and it might differ in phase by a certain angle (p a - (For example: wave re- 
flecte by the ffee-end of the string is known to have a phase change of n.) The real amplitude 
A and the phase angle f can be combined into a complex amplitude C = Ae'^ . 

Note that — x direction can be viewed as performing a coordinate transformation of the xy 
coordinate system by n; that is, 

Q—ikx __ gikr cos(0+7r) (82) 

Then, this reflecte planar wave can be approximated as, according to eqn. (73), 

OO OO 

Ce~ ikx = Ce tkrcos ^ +n) = C ^ i n J n (kr)Q in ^ d+n) = C ^ (-i) n J n (kr)e in0 (83) 

n = — oo n=—oo 

For individual scatterers, the wave expansion coefficient for the propagation component 
of the reflecte wave can be expressed as 

{Af}„ = Ce~‘ kx '(—i) n (84) 


V.4 Wave Expansion for a Point Load 

When a concentrated force F of wave number k is applied at a location denoted as x 0 , the 
resulting displacement is expressible as (L. Horchens, Imaging of Materials Inhomogeneities 
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with Flexural Waves, Ph. D. Thesis, Delft University of Technology, The Netherlands, 2010, 
Chapter 4) 

WinC(x) = -? 8 (k lX ~ * o|) " H { ~ lk |X “ x °^ (85) 
Note that in the derivation for the above expression (in Appendix B), the time-incorporated 
expression for a general forward-traveling wave is fit — x/c ) for one-dimensional case. 
Consequently, a harmonic wave would be as Q l (a>t-kx) ^ w ^ ose time factor is e l0)t . In the 
two-dimensional case, the spatial factor for the an outgoing cylindrical wave in a polar co- 
ordinate system would be Hf'(-), as the asymptotic expressions for the Hankel functions 
H n),<2> (x) = Ce ±lx where C is a coefficient the positive sign is for the Hankel function of 
the firs kind, and negative for the second kind. In other words, the Hankel functions of the 
second kind in eqn. (85) represent out-going waves. This notation differs from most contem- 
porary literature, in which Hf\-) is generally used to represent an outgoing wave, whose time 
factor, if to be included, is e~ lcot . In adopting the above expression for the present analysis, 
Hf\-y s are simply replaced by the corresponding Hff-y s with the same arguments. 

Furthermore, if a local coordinate system ( r s , 9 S ) is constructed such that its origin O s is 
located at the point of incident wave, the above expression can be simplifie to 

= ~'Wd “ 

= ~WD [ tH »" (kr ^ + \ K °( kr f < 86 > 

where the relation between Hf\-) and K 0 (-) in eqn. (11) have been used. Note that the inci- 
dent fiel is omnidirectional, and is thus independent of 6 S . 

In most computations for scattering problems, all the wave amplitudes are normalized 
by that of the incident wave. In such cases, the constant coefficien in W mc (rj, 9j) becomes 
immaterial and can thus be dropped. On the other hand, since the response due to a concen- 
trated force is a singular field the normalization factor can be the fiel amplitude at a certain 
characteristic length a from the source, as 


where 


W (r s ) = 


W™(r s ) 

|IU inc (tf)| 


1 

A 


iHq \kr s ) + — K 0 (kr s ) 

71 


(87) 


A = 


iHi fka) H Ko(ka ) 

7 X 


( 88 ) 


For convenience and in keeping consistency with the case of planar incident wave, the nor- 
malized incident wave is simply taken as the incident wave, and the overbar in the above 
expression will be dropped in the discussions to follow. Using a wave expansion basis based 
on this source-local coordinate system, the matrix notation for the point source is 

Wr(r s , 9 S ) = {A?} T {H{r s , 9 S )} + { A e s } T {K{r s , 9 S )} (89) 
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Then, coefficient are such that 

i 2 

{Af } o = -^ {v4f}o = -— {Af}„ = 0 for other n (90) 

In scattering computations, it is necessary to express the incident wave in individual scat- 
tered’ local coordinate system, say (r, , 6j). Then, coordinate transformations wave wave 
expansion bases in eqns. (58) and (69), 

{A') T {H(r„e,)} = {A!) r [R p s ,]{nr u e,)) (91) 

{A’} T {K(r s , «,)} = {Ai) T \g^\{I(r t , «,)} (92) 

where [Xfi] and | R’i ] are coordinate translation matrices for the propagating and evanescent 

components between local coordinate systems ( r s ,6 s ) and (r,-, 6 A, respectively. Note that 
there is only one non-zero entries in the coefficien matrices { A f } and {A*}. The matrix 
multiplication can be carried out as zero-th column of the respective [ R ] matrices. If the point 
force incident wave is expressed as 

w; nc (r t , e t ) = {Af} T {J(n, Oi)} + {AlfiHru 9,)} (93) 

the wave expansion coefficient for this incident wave are 

{Af} n = ~e~ ine °‘H- n (kd si ) = (-1 ) n+l ^e- ind °<H n (kd si ) (94) 

{A e i)n = -2-(-l ) n Q -^K_ n (kd si ) = — ?-(-l) n <r tn9 ' l K n (kd S i) (95) 

7T A n A 

where (d si , 6 S i ) are the coordinates of the source located in the (r,-, 9i) coordinate system, and 
the relation K- n {z) = K n (:) has been used. 

Note that if the evanescent component can be neglected, the propagating component of the 
incident wave is the same as a point source in a two-dimensional free space. 

In the case of multiple wave sources, then only one of the wave will be used to defin 
the normalization factor as described above, and all remaining sources can have a relative 
amplitude and and a phase difference from the one chosen for the normalization. The phase 
difference is simply expressed by e'^. 

V.5 Wave Expansion for a Line Source 

A line source can be constructed from an array of identical point sources lined up along the 
designated segment. This makes sense both physically and mathematically. Mathematically, 
the sum becomes a Schlomilch series if the line is infinitel long, which has been shown 
by Twersky to turn into (expressible as) planar waves for the propagating components. The 
complication comes from the normalization, especially when the length of line is finite This 
will be dealt with numerically. In the numerical computations, The wave fiel due to the 
multiple sources will be calculated along a line a distance a from the line source. The average 
amplitude of the incident wave along the line segment will be used as the normalization factor. 
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VI. Computational Considerations 


■ Since the Young’s modulus E does not enter in the formulation except in computing the 
rigidity D, which have different expressions for a uniform plate and a sandwich panel, in 
the code, E is not a required material property. Instead, D is used as an input. D can be 
calculated using a spread-sheet external to the code. 

■ The thickness h and mass density p are needed in relating the wave frequency to the wave 
number. Note that when the plate is a composite panel, h is the total thickness, and p must 
be the overall mass density. In other words, ph should forms the areal mass density of the 
composite panel. 

■ Poisson’s ratio is extensively used. Thus it is needed as an input. 

■ Since the wave speed is frequency-dependent, the wave number no longer represents a sim- 
ple (linear) normalization of the frequency. Since it is desirable to have the computational 
results presented on a linear scale of the frequency, the wave number k is not used as a 
normalized frequency in the computations. For the lack of a better alternative, the Hertz 
frequency / = co/(2tz) is used as the input. Correspondingly, all other parameters will be 
un-normalized, specifie in SI units. 

A sensible alternative for normalizing the frequency would be to define 


such that 




co 
co o 


= (kh) 2 


(96) 


(97) 


This way, co/coo = 1 means kh = 1, and in turn the wave length is 2 n times of the plate 
thickness. 

■ In multiple scattering formation, only [/? /J ] and [R e ] are used (that is, not the regular counter 
part) for coordinate transformation. They have a unifie expression expressible as 

[Rij]nm = (-\) m Q l(n ~ m)e U k? n _ m (kd,j) (98) 


where 3f„(z) will be H n (z ) for [R p \, and K n (z) for [R e ]. In other words, they can be 
computed using a unifie function. 

■ The following product of generalized matrices are used in the multiple scattering analysis: 


MMy] 7 


- [T P P] 

[7T1 ■ 

1 

C: 

i — 

o 


■ [TP ][K£] 

1 

(j ' 

. [T- P ] 

IV 1 . 

i — 

o 

[*,q j 


. [T‘ P V„] 

1 


(99) 


It is clear that directly constructing the generalized based on the right-most expression is 
more efficien than constructing the two matrices separately then performing the product of 
the generalized matrices. 
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■ The wave fiel can be obtained directly by the product of the generalized matrices for the 
wave expansion coefficient and bases. Thus the multiple scattering solution, that is, the 
wave expansion coefficient for individual scatterers, will be stored in the generalized ma- 
trix form. 

■ For an elastic scatterer, it is assumed that the thickness is the same as the host plate. Uneven 
thickness will causes three-dimensional effects and would be beyond the capability of the 
current model. 


VII. Validations 

VI. 1 Energy Flux Conservation Validations 

According to Norris & Vemula (Scattering of Flexural Waves on Thin Plates, JSV, 118 ( 1 ), 
115-125, 1995), the energy flu conservation for the fl xural wave is expressible as 


coD 


-3 


3v 2 ff a w 
w — - v 2 w- 


' c 


dr 


dr 


ds > = 0 


( 100 ) 


The physical meaning of this relation is the principle of energy conservation: if the surface C 
(a path in a two-dimensional problem) enclosed neither a source nor a sink, the energy enter 
into the surface should equal to the energy emit out of the surface. 

In general, the wave fiel on the plate can be divide into propagating and evanescent com- 
ponents, 


w = w p + w e 

and they satisfy different equations: 

W 2 W P + k 2 W p = 0 


and 


W 2 W e -k 2 W e = 0 


or 


V 2 W P = —k 2 W p and 
The integral in eqn. (100) can be expressed as 


V z W e = k 2 W e 


( 101 ) 

( 102 ) 

(103) 


-L 

=k2 L 

=k2 l(~ 


dV 2 W 


- V 2 W 


rdW 

dr 


ds 


(W p + W e ) 


d(-WP + W e ) 

dr 


- {-WP + W e ) 


d{W p + W e ) 


dr 


ds 


JW dW dW dW 

w p — — + w p — w e —~ + w e 


dr 


dr 


dr 


dr 


— p dW p — P dW e — e dW p — e dW' 


+ W 


dr 


+ W 


dr 


W 


dr 


W 


dr 


) 


ds 


(104) 
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(105) 


And in the far field only the propagating component remains, and the inegral becomes 

f ( dW P — p dW p \ , 

/farfiel = k J l~ W + W ~^~) dS 

Note that the above relation must hold for all incident waves. It is possible to construct an 
incident wave that constant only one entry in the wave expansion, leading to a pair of scattered 
waves. For example, if the incident is An = 1, the resulting total wave fiel is 

W p = J n (kr) + T pp H n (k r)e in0 W e = T n pe K n (k r)e in6 (106) 

If the integration path is circular with a large radius (far-field) eqn. (100) gives 

0 = 3 {(J„(kr) + T™H n (kr )) + TZK(kr)) 

- (j n (kr) + T P n P JZ(kr)) (J' n (kr) + T"H'„(kr))\ 

= 3 { J„(kr) [TZK(kr) - T™H' n (kr)\ 

+K(kr) [T™ H n (kr) -TZjZ(kr)\ 


+ ra 2 [H„(kr)H^kr) - H'„(kr)H n (kr)\ j (107) 

Note that, since J n (kr) and J'nQcr ) are real, 

3 \jZK(kr) - T™H'„ikr)\ = -23 {T™H' n (kr)} (108) 

3 \l™H n (kr) - TZjZikr )] = 23 {T™H„(kr)} (109) 

The equation becomes 

3 {277/ [ J„(kr)H'„(kr ) - J'„(kr)H„(kr )] 

+|7’/f| 2 [H„(kr)H a >' n (kr) - H' n (kr)H™ (*r)]} = 0 (110) 

Recall the Wronkian relations for Bessel and Hankel functions in eqn. (31) and 

H£\z)H m ' n {z) - H“>'(z)H?>(z) = — (111) 

71 Z 

Then, 

3 { 2T "^ + ' T "^} = 0 (112> 

which can be simplifie to 

9t{r//} + |r//| 2 = o (ii3) 


In the near field noting K n (kr ) is real, 


dW —edw e 
W e — W — = 0 


dr 


dr 
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I = 2 nk z {TZK(kr) [J,(kr) + T™H„(kr)] - T%K„(kr) [j' n (kr) + TZX(kr) 
+T„ZK' n (kr ) [j n (kr) + r"7/„(A:r)] - T P n ‘ n KAkr) [j' n (kr) + TffH^kr)]} 

= T£ [ K'„(kr)J n (kr ) - K n (kr)J'„(krj\ + T% [ K'„(kr)J n (kr ) - K„(kr)J'„(krj\ 
+T£TZ [ K' n (kr)H„(kr ) - K„(kr)lf n (kr)] 

+TZ„T [ K' n (kr)H„(kr ) - K„(kr)H' n (kr)} 

= 2» {r//} - K n (kr)J' n (kr)} 

+23t [t^tZ [K'„(kr)H„(kr) - K n (kr)H'„(kr )]J (114) 

and the condition is automatically satisfie without posing any conditions on the T -matrices. 
Similarly, if the incident is A e n = 1 , the resulting fiel is 

W p = T^H n (kr)e ind W e = I n (kr ) + T“ K n (k r)e ind (115) 

In the far fiel with only the propagating component, 

[ r ( jWp — 3 w p \ , ) 

o=3 {icr ^ ) 

= 3 {T;!H n (kr)r n P X(kr)-rZH„(kr)T;i;H:(kr)} 

= 3{|rn 2 J 7 } (116) 

The content in the pure imaginary. This would require | T P n \ = 0, which is not correct. This 
suggests that the condition is not satisfied This is likely due to the fact the the incident wave is 
not included; meaning that the path actually include the source. The source is the evanescent 
incident, which is localized in the vicinity of the origin. 

Thus, to validate the other part of the calculation for the T -matrices, it would be necessary 
to numerically compute an integral over a closed path that avoids the origin. 

VI.2 Optical Theorem Validation 

The scattering form factor is the far- fiel angular distribution pattern of the scattered wave 
when it is subjected to an planar incident wave of unit amplitude. It is define as (Norris & 
Vemula, JSV 1995) 

lim IV scr = _L eKkr-n/A) f^ (117 ) 

r^oo V2 r 

For the single scatterer problem, recall the asymptotes for Hankel functions and modifie 
Bessel functions 

H " (kr) = 


J(kr—nn/2—n/4) __ ^—r 


j(kr—n/4) 


( 118 ) 
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and 


which means that the evanescent weaves do not contribute to the far field Thus, eqn. (17) 


OO I Z OO 

lim W SCI = V BEH„(kr)e ind = J _^_ e f( * r-,r/4) V B% 

OO n n\ , y ^ n 


i~ n e n 9 


( 120 ) 


m = 7^1 

The scattering cross section can be calculated as (Norris & Vemula, 1995) 


( 121 ) 


°=2j 0 \WW de 022) 

In the mean time, the optical theorem for the fl xural wave is (Norris & Vemula, 1995) 


(123) 




which is the consequence of the energy flu conservation in the far-field Thus, combining the 
above two equations gives the following identity that is equivelent to the energy flu conser- 


vation: 


For a single scatterer problem, the integration in eqn. (122) gives 


(124) 


<7 = 1 O 1 " f(6)f(6)de = ~ 


C 2n / oo \ / oo 
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\n = — oo / \m=—o o 


Zi m Q- im0 do 


r 2n oo oo 

1 EE 
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B?Bti m - n Q i ( n ~ m)e d6 
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r\ OO , OO 

- E = r E KT 


(125) 


where when expanding the two summations, noting that 


i(n—m)Q dQ — 


2n when m = n 

0 when m ^ n 


(126) 
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(127) 


On the other hand, the evaluation of eqn. (123) gives 


a = — 2 



Comparing the above results we have 


4 

k 


E 


OO 

E (|^| 2 + 9i{^r«}) = 0 (128) 

n——oo 

For the multiple scattering problem, the optical theorem is still valid. Thus, the equivalent 
validation in eqn. (124) remains valid. It boils down to computing the scattering form factor. 
One approach is to perform a similar “scatterer polymerization” scheme to turn the scatterer 
ensemble into one abstract scatterer. Alternatively, it can be performed numerically. In the 
numerical computations, the evanescent mode can be neglected. 


Appendix A. Definitions for ©-Functions 


A.l For Propagating Waves 

Consider a serial expression for a deflectio is expressed as 

OO 

W(r) = J2 oi n ^n{kr)e nd 

n——oo 

where % t (") the either J„(z) or H n (z), and z = kr will be used as an abbreviation. 
Recall that in a polar coordinate system, 

d 2 W 13 W 13 2 W 

dr 2 r dr r 2 3 6 2 


and 


Then, 


i dw 

r dr 
1 d 2 W 

V 2 ~d9 2 ' 


J OO OO 

- u n tf'(kr)e tnd = z J2 °^ r K(kr)e in9 

n=—o o n = — oo 

OO 

- E 

n = — oo 


d 2 w 


2 


VlW = + E 7 T [~K(z) ~ n 2 V n (z)\ e 


JnO 


(A.l) 


(A.2) 


(A.3) 

(A.4) 


(A. 5) 
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On the other hand, W satisfie the Hel mh oltz equation, that is 

(X) (X) 

V 2 W = —k 2 W = — k 2 ^ a n c to n {kr)Q ine = — z 2 ^ ^■% t (z)e i 

n=—oo n = — oo 

Combined, they give 


in 6 


(A. 6) 


00 ry 

^r= E ^ [-^W + (» 2 -**)%(*)] e M 


(A.7) 


Thus, following the definition for M rj Me and K,., 

d 2 w' 


Mr = ~D 


vV z W + (1 - v) 


dr 2 


JnO 


= ~D [-vzV„(z) - (1 - v)z^(z) + (1 - v)(n 2 - z 2 )tf n (z)] e l 

n——oo 

00 

= ~D ^ ^ {-(1 - v)z c tf' n (z) + [(1 - v)n 2 - z 2 ] tfnikr)} e ine (A.8) 


M e = -D 


V 2 W -( 1 -v) 


d 2 W' 

dr 2 


-D ^ -^[-z 2 tf n (z) + (l-v)z^(z)-(l-v)(n 2 -z 2 )tf n (z)]e ine 
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00 
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(A. 10) 

(A. 11) 
(A. 12) 
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Then, 


V r = ^2 a„ [-z 3< ^(z) - (1 - v)(—n 2 ) < rf„(z) + (1 - v)(-n 2 z)%(z)] e in 9 

n——oo 
D 00 

= -— ^ a„ {- [(1 - v)n 2 + z 2 ]z^(z) + (1 - v)n 2< rf„(z)} e‘ nd (A.13) 

n=—oo 

A.2 For Evanescent Waves 

Consider a serial expression for an evanescent wave as 

OO 

W(r) = a n £?? n (kr)e ind (A. 14) 

n =— oo 

where 3f n (z) is either I n (z) or K n {z). The evanescent wave satisfie the modifie Helmholtz 
equation, as 


V 2 W = k 2 W = z 2 °^^ n (kr)e 


inO 


Thus, other components are as the following 

If = £ £ 


inO 
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n=—o o 
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n — — oo 
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dd 2 

d 2 w 
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A (V 2 ^) = z J £ ^„'(z)e« 


n=—o o 
oo 


a 3 w 
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-n 2 - a„^; / (z)e ? 




Thus, following the definition for M r , Mg and K,., 

a 2 w 


Mr = -D 


vV 2 W + (1 - v) 




(A. 15) 


(A. 16) 
(A. 17) 
(A. 18) 
(A. 19) 
(A.20) 


JnO 


= —D ^ [vz 2 J^(z) - (1 - v)zJ^'(z) + (1 - v)(n 2 + z 2 )J^(z)] e z ‘ 

n——oo 

oo 

= -D {-(1 - v)z^'(z) + [(1 - v)n 2 + z 2 ] 3? n (kr)} e ,n9 (A.21) 
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(A.22) 



In a unifie notation, using 3F* (•) to represent any of Bessel functions. Specificall , 

&nQ = U 0 ^ 2 Q = H n (.) = /„(.) 2f n \-) = K n (.) 

If a displacement is expressible as 

OO 

W(r) = oi n ^(kr)e ne (A.24) 

n = — oo 

Then, the resulting moments and shear force can be expressed by a series of ©-functions such 
as 

| OO 

M r = — - ^ a„©j (n,kr)Q ine (A.25) 

n — — oo 
I oo 

Mg = — - ^ a„©2 (n,kr)e in9 (A.26) 

Yl — — OO 
1 OO 

ff = — 3 ^2 oi n & 3 (n,kr)Q ind (A.27) 

n——oo 

These ©-functions are define as 

©*,(«, z) = D j-(l - v)zJ^'(z) + [(1 - v)n 2 =F z 2 ] 2%(z ) J (A.28) 

& 2 (n,z) = D j(l - v)zJ^'(z) - [(1 - v)n 2 ± vz 2 ]J^(z)j (A.29) 

& 3 (n,z) = D j- [(1 - v)n 2 ± z 2 ]z2?*' n (z) + (1 - v)n 2 J^(z)} (A.30) 

where the upper sign in front of z 2 is taken when t = 1 or 2, and the lower sign is taken when 
t = 3 or 4. 
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